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INFLUENCE OF MEAN FLOW ON 
BOUNDARY LAYER GENERATED INTERIOR NOISE 
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SUMMARY 

An analysis is made of the "interior noise" produced by high, subsonic 
turbulent flow over a thin elastic plate partitioned into "panels’* by straight 
edges transverse to the mean flow direction. This configuration models a 
section of an aircraft fuselage that may be regarded as locally flat. The 
analytical problem can be solved in closed form to represent the acoustic 
radiation in terms of prescribed turbulent boundary layer pressure 
fluctuations. Two cases are considered: (i) the production of sound at an 

isolated panel edge (i.e., in the approximation in which the correlation 
between sound and vibrations generated at neighboring edges is neglected) , and 
(ii) the sound generated by a periodic arrangement of identical panels. The 
latter problem is amenable to exact analytical treatment provided the panel 
edge conditions are the same for all panels. Detailed predictions of the 
interior noise depend on a knowledge of the turbulent boundary layer wall 
pressure spectrum, and are given here in terms of an empirical spectrum 
proposed by Laganelli and Wolfe. It is expected that these analytical 
representations of the sound generated by simplified models of fluid- structure 
interactions can used to validate more general numerical schemes. 
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1. INTRODUCTION 

Interior cabin noise of passenger aircraft is usually attributed to 
airborne sources , which include the direct incidence of engine generated sound 
on cabin sidewalls and the high speed turbulent air flow over the exterior 
fuselage, and structure borne sound (transmitted to the cabin via vibrational 
motions of the airframe) caused by engine vibration, and by vibrations induced 
by impingement of engine exhausts on wings and flow control surfaces. Flight 
tests [1] with engines "off" indicate that the boundary layer is responsible 
for a substantial component of the interior noise. It is significant for 
practically all types of aircraft, and is usually considered to be one of the 
most important sources of cabin noise for a jet powered passenger aircraft in 
steady flight. The spectral levels and the effective frequency range of 
boundary layer pressure fluctuations change with increasing subsonic flight 
speed [2], and boundary layer generated cabin noise can vary significantly 
with frequency, in contrast to the smoother variations of the exterior wall 
pressures [3]. This is because fuselage transmission characteristics and 
boundary layer scattering mechanisms are critically dependent on the geometry 
of structural changes in the cabin walls. 

To control interior noise one must be able to predict both the amplitude 
and phase of the cabin sound pressure field. For the random noise sources in 
a turbulent boundary layer, the most that can normally be done is to determine 
the space- time correlation of the interior sound, and possibly also the 
correlation between vibrating structural elements and the sound. Turbulence 
generated sound enters the cabin both by direct transmission through the 
fuselage and via the excitation of structural vibrations that subsequently 
generate sound. A statistical analysis of structural excitation and noise 
transmission must include the influence of the space- time correlation of the 
boundary layer wall pressures. Most work so far conducted on this problem has 
taken little or no account of the details of the spatial characteristics of 
the wall pressures [4] , although it is well known that the convection of 
strong "hydrodynamic" pressure fluctuations over the fuselage is important 
when "coincidence" occurs (when the convection velocity is equal to the phase 
speed of a resonant wall mode) , causing a large increase in the interior noise 
and wall vibrations [5, 6]. 
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Hitherto it has been usual to simplify the analysis of such problems by 
neglecting the coupling between the structural and acoustic equations of 
motion. Relatively little work has been done on the fully coupled equations 
[7] in aeronautical applications, although the validity of the decoupled 
approximation is suspect at higher subsonic Mach numbers, when dynamic 
boundary layer pressure fluctuations can be large, and it is known to fail in 
analogous underwater applications where the coupling is large because of the 
relatively high density of water [8]. 

In this chapter an analytical investigation is made of a simplified model 
of the generation of interior noise by subsonic, high speed turbulent flow, 
when the fuselage is assumed to be locally flat. The coupled, fluid- structure 
equations can be solved in closed form in two general cases. In the first of 
these we examine the sound produced by turbulent flow interacting with a 
single straight edge, transverse to the mean flow, that separates two large 
elastic panels. This is a situation in which the sound and vibration produced 
at neighboring edges can be regarded as uncorrelated; it is unreasonable in 
practice, because it neglects strong couplings between neighboring edges via 
flexural motions of the panels. Next, this problem is extended to the case of 
periodic, identical panels with transverse edges, whose interaction with the 
boundary layer can be solved in closed form provided the panel edge conditions 
are the same for all panels. This arrangement may be regarded as a local 
model of periodic fuselage ring stiffeners, which are usually the dominant 
restraints on panel motions. Longitudinal "stringers" between the stiffeners 
are important in controlling panel resonance frequencies, but, because they 
are parallel to the exterior mean boundary layer flow, their contributions to 
the generation of sound and vibration by fluid-structure interactions will 
tend to be of less significance. Similarly, recent acoustic -holographic 
measurements [9] indicate that the more visually obvious structural 
discontinuities associated with cabin windows are also unlikely to be 
important sources of noise and vibration. 

The analytical model is formulated in §2, where the reverse flow 
reciprocal theorem (discussed briefly in the Appendix) is applied to determine 
a Green’s function for treating the "scattering" of boundary layer wall 
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pressures at panel edges. This is used (in §3) to investigate scattering at 
an isolated straight edge, when the edge is regarded as either clamped or 
simply supported . The edge radiated sound is expressed in terms of the 
wavenumber -frequency spectrum of the boundary layer wall pressure fluctuations 
[10, 11], whose properties are specified independently. Here we make use of 
an empirical model proposed by Laganelli and Wolfe [12]. The principal 
characteristics of the sound produced for either edge condition are similar, 
being dominated by the requirement that the panel displacement f , say, must 
vanish at the edge, rather than conditions related to restrictions on higher 
derivatives of f . The more general problem of a periodic array of coupled 
panels is discussed in §4. As the frequency varies, the coupling between 
neighboring panels causes large fluctuations in sound levels about those 
predicted for the isolated edge. The analysis is expected to be valid when 
the Strouhal number co<5*/U < 1 (w being the radian frequency, 6* the boundary 
layer displacement thickness, and U the velocity of the main stream), and we 
shall actually conclude that, for subsonic mean stream velocities, it is only 
when this condition is satisfied that significant production of sound occurs 
by the fluid- structure interaction. The results should be of value in 
validating more general numerical schemes for interior noise prediction. 
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2. THE GOVERNING EQUATIONS 


2.1 The aerodynamic sound problem 


Consider high speed turbulent flow of fluid of mean density p l and sound 
speed over the surface of a thin elastic plate that occupies the plane x 2 “ 
0 of the rectangular coordinate system (x^x^Xj) (see Figure 1). The mean 
flow is in the positive x : -direction in the "exterior" region x £ < 0 . The 
plate is partitioned into panels by rectilinear edges at x x — X x , X 2 , 

X 3 It is required to determine the sound radiated into the "interior" 

region x £ > 0 as a result of the fluid- structure interaction. The fluid in x 2 
> 0 is at rest in the undisturbed state, and has mean density and sound speed 
respectively equal to p Q and c Q . 


Take the equations of aerodynamic sound in the form [13] 


1 3 2 

c z dt z 

1 fd a 3 2 
c^lat + u ax 1 J 


V 2 j B — 0 , x 2 > °- 

V 2 ]b - S(x,t), x 2 < 0. 


( 2 . 1 ) 


( 2 . 2 ) 


In these equations B = w + ^v 2 is the total enthalpy, w being the specific 
enthalpy of the fluid, v the velocity, U is the uniform main stream velocity 
in the exterior region, and t denotes time. The aeroacoustic source term 
S (x , t) vanishes except in those regions where the vorticity and entropy 
gradients are non- zero, but its functional form will not be required in the 
following analysis. When B is known, the acoustic pressure p can be 
calculated from the formula 


I (2.3) 

p at Dt 

where p is the local mean density and D/Dt is the material derivative. This 
equation is valid at high Reynolds number provided the motion may be regarded 
as isentropic (i.e., the entropy of each fluid particle is constant), which we 
shall assume to be the case. Equation (2.2) is applicable provided the 
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characteristic wavelength of the sound is large compared to the boundary layer 
displacement thickness 6* , which should be the case at radian frequencies a> 
satisfying /\J < 1. 


The panel motions are assumed to be governed by the linear thin plate 
equation, such that if f^.x^t) is the flexural displacement in the 
x 2 -direction , then 

(dV* + m^Jf + [p] = 0, (2.4) 

where V* = (d z /3x 2 + d 2 /dxl) 2 , D is the bending stiffness of the plate, m the 

2 13 

mass per unit area, and 

[p] = P(x x ,+0,x 3 ,t) - p(x x , -0,x 3 ,t) (2.5) 

is the pressure loading. The linearized approximation is expected to be 
adequate for aircraft panels, where displacement amplitudes do not usually 
exceed 5 ^m. In practice panel motions are damped by interior "trimming" to 
reduce interior noise levels. We shall take approximate account of this by 
introducing a panel loss factor r ) , such that for vibrational motions of radian 
frequency u> (proportional to e i<Jt ) , 


D - D o [l - i»7 (cl>) ] , r7(-w) = -»/(w), 


( 2 . 6 ) 


where D > 0, and n > 0 for w > 0 [14]. The panel and fluid motions are also 

O 

coupled by the linear relations 


d z £ 

at* 


1 dp 


x 2 = +0; 



1 dp 


X„ 


- 0 , 


(2.7) 


which follow from the linearized, normal components of the fluid momentum 
equation just above and below the plate. 


2.2 Green’s function 

Let G(x,y, t-r ) denote the time domain Green’s function defined by the 
solution of equations (2.1), (2.2) with outgoing wave behavior when the source 
term S in (2.2) is replaced by an impulsive point source 6(x-y)6(t-r) at y. 

The acoustic pressure p(x,t) in the interior region (where B » P/P 0 ) is tlien 
given by 

p (x , t ) = J G(x,y,t-r)S(y,T )d 3 ydr , (2.8) 
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where the spatial integration is over the boundary layer sources and the time 
integration is over -® < r < Green’s function can be partitioned in the 

following way: 

G(x,y , t-r ) =G o (x,y,t-r) + G s (x,y , t-r ) , (2.9) 

where G is Green’s function for a uniform, homogeneous plate, and G is the 

O S 

additional component attributable to the panel edges. 

In practice the acoustic pressure p(x,t) is required in the interior 
domain at large distances from the source region, typically at distances |x-y| 
from the sources that are large relative to the acoustic wavelength. This 
observation permits the derivation of a simplified form of the Green’s 
function . 

To do this set 


G( 


x,y , t-r ) = -if G(x,y,w)e' iw ^ t ' r ^ 
2?rJ _ 


do) , 


( 2 . 10 ) 


where G(x,y,ci;) satisfies 


[V 2 + /c 2 ]G(x,y,w) =0, x 2 > 0, 

[V 2 + (k 1 + iM3/3x 1 ) 2 ]G(x,y,w) = 6(x-y), x 2 < 0, 

and 

«o - “/ C o . “ w / c ! . M “ U / C ! • 


( 2 . 11 ) 


( 2 . 12 ) 


For observer positions x many acoustic wavelengths 2x//c 0 from the point 
source at y, we invoke the reverse flow reciprocal theorem, according to which 


G(x,y,u) =— G R (y,x,w), (2.13) 

'o 

where G R (y,x,w) is the solution of equations (2.11) when x and y are 
interchanged (so that the source is now regarded as being at x, in the 
interior domain, and the observer at y in the exterior flow region) and M is 
replaced by -M (i.e., the mean flow velocity is reversed; see Figure 2). This 
result is a particular case of a more general reverse flow reciprocal theorem 
(see, e.g., [15, 16]) whose proof is outlined for the present problem in the 
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Appendix . 


When | x | -* ® the sound wave G R (y,x,u>), say, generated by the reciprocal 
source at x may be regarded as a plane wave impinging on the elastic plate. 
Thus, to calculate the reciprocal Green’s function G R in this limit we have 
only to solve a plane wave diffraction problem, in contrast to the direct 
calculation of G, which involves the interaction of the near field of a point 
source at y with the inhomogeneous plate. The incident, reciprocal plane wave 
(regarded as a function of y) is evidently given by 


Gj (y,x, w) 


-e^o i x *yi 

47r|x-y| 


- e iK o { l x l ~ n *y } , |x| - CO, (2.14) 

47T | x| 


where n is the unit vector 


n = x/|x| . 


(2.15) 


The diffraction problem is solved in two stages. First we determine the 
motion on and near the plate produced by the incident wave when the panel 
edges are ignored, i.e., for a homogeneous plate. This defines the reciprocal 
Green’s function G R (y,x,w), say, that determines G o (x,y,t-T) of (2.9). 
Recalling that G R (y,x,w) represents the perturbation total enthalpy, it is 
easy to show in the usual way (by making use of the homogeneous forms of 
equations (2.11) in the reciprocal problem and of equations (2.3), (2.4), 

(2.6) [17]) that the displacement f of the homogeneous plate produced by the 
reciprocal source can be expressed in the form 

(2.16) 


where 


3(n,w)e 1,c o l x l 
2»riuc 0 | x | 


(2.17) 


3(n,w) 


ie 




( n i +n 3 




(PoCo/PiCi)(c 0 /c 1 -Mn 1 ) 2 ~|3 

V{ ( c 0 /<=i -Mn 1 ) 2 - (n 2 +n 2 }J J 


(2.18) 


In these expressions 

V{(c 0 /c 1 -Mn 1 ) 2 -(n2 + n^)) _ sgn ( Co/Ci ) | (c 0 / Cl -Mr^ ) 2 - (n z +n z ) | 1/2 

= isgn(w) | (c 0 /c 1 -Mn 1 ) z -(n 2 +n 2 ) | 1/z 
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according as the argument of the square root is < 0. The quantity u> c (> 0) is 
the coincidence frequency c z (m/D o ) 1/z , at which the in vacuo phase velocity of 
bending waves on the plate in the absence of damping is equal to the speed of 
sound in the interior region. Then e is a fluid loading parameter defined by 

€ _ P ° C °_ . (2.19) 

u) m 


For an aluminum plate in air, for example, e » 0.0022, and w.h/c 0 ** 0.22, 
where h is the plate thickness. 


To complete the diffraction calculation it is necessary to examine the 
interaction between the panel edges and the incident flexural disturbance 
(2.16). This interaction produces scattered flexural motions on the panels, 
which we denote by f s (y : ) e' ln 3 y 3 , and potential flow motions in the fluid 
which satisfy the reciprocal forms of equations ( 2 . 11 ) (with the source 
omitted). If the velocity potentials are denoted by <p Q e in 3^3 and cp 1 e in 3^3 
respectively for y 2 < 0 , it follows from the reciprocal equations of motion 
that , 


CD 

f s (y x ) = J f s (k)e lk ^i dk, 
- 00 



( 2 . 20 ) 


where 7 = y(k) - [/c 2 ( 1 -n 2 ) -k 2 ] 1/2 (branch cuts being chosen such that, for real 
k, sgn(y) - sgn (u>) when 7 is real and is otherwise positive imaginary). The 
function T + (k) and a related function T_ (k) , also needed below, are defined 
for real k by 

T ± (k) = sgn(/c 1 ±Mk) | (/^iMk) 2 -/c 2 n 2 -k 2 | 1/2 , for (»c 1 ±Mk) 2 -«; 2 n 2 -k 2 > 0, 

= i | (/c 1 ±Mk) 2 -/c 2 n 2 -k 2 | 1/z , for (/Cj ±Mk) 2 -/c 2 n 2 -k 2 < 0 . ( 2 . 21 ) 
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Using the representations (2.20) and (2.3) (noting that B — iuxpj • j = 0, 
1 in irrotational flow), we find in the reciprocal problem that the jump in 
the scattered pressure p s e' 1 ' c o n 3 y 3 across the plate is given by 



2 

P o w 
7(k) 


py (a>+Uk) 2 1 (k) iky 

T + (k) P*' 


dk, 


( 2 . 22 ) 


where [p ] is defined as in (2.5). The scattered pressure and displacement 
must satisfy the thin plate equation (2.4) except at the panel edges, i.e., 
substituting from (2.20) and (2.22) into (2.4), we must have 

00 

J jr + (k, w)f s (k)e iky i dk = 0 , y 1 * X 1 , X 2 , X 3 (2.23) 


where 


<e ± (k,u) 


D ( »c 2 n 2 +k 2 ) 2 - mw 2 


W(k) 


Pi (wtUk) 2 "! 

r ± (k) J ’ 


(2.24) 


Now jC + (k,u>) ■ 0(k 4 ) as 
finite as y 1 •» ^±0 (N = 1 , 
taken in the form 

r, (k> - E„ r, 


| k | -» oo, which implies that if f (y ) is to remain 
2, 3 ,...), the general solution of (2.23) must be 


[\o +A N1 k+ V k2+A N3 k3 ) 


e’ ikX N 
JC + (k ,w)’ 


(2.25) 


where the coefficients A^ = 
determined from the panel edge conditions. Thus, if f - Z(y x )e 1 ^ c o n 3^ r 3 is the 
net displacement of the plate in the reciprocal problem, including the 
component f of (2.15) for the homogeneous plate, we can write 


(u>) are independent of k, and are to be 






dk + e‘ i/c o n iyi. (2.26) 


It remains to determine the values of the coefficients from 

conditions specified at the panel edges y 1 « X 1 , X 2 , X 3 The details of 

this calculation are outlined in §§3, 4 for two different panel 
configurations. The velocity potential tp 1 defined in (2.20) may then be 
evaluated by substituting for f (k) from (2.25). Then (2.10) and the 
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reciprocal formula (2.13) yield the following representation for the component 
G of the time domain Green’s function 


G s (x,y, t-r) 


2 *P 0 


j"« Vl (y lf «)e- i I ,6 o n 3y3 +w ( t - T >] dc, 

- 00 

y 2 < o, |x| 


00 , 


(2.27) 


Here the source point y (y < 0) is located within the mean flow in the 
exterior domain, and the observer position x is in x £ > 0 in the interior 
region. Substituting for < -p we find 


G (x.y.t-r) 


-P l 


( 2rr) 2 c Q p 


J f + A k+A 2 k> + A_ kd 

0 1 x | r + (k)r + (k lU )P 0 ^ j 


x e i{k( yi- x N ) - r + y2-' c o n 3y3' a,(t ' r 'i x i /c o )> dkd W , 

y 2 < 0, | x | -» <». (2.28) 

The condition that G s is a real valued quantity implies that 

Agj (-«) = (-D^j <«). (2.29) 

where the asterisk denotes the complex conjugate. 
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3. NOISE GENERATED AT AN ISOLATED PANEL EDGE 


3.1 The acoustic pressure spectrum 


To assess the importance of different edge conditions we first assume the 
infinite plate of Figure 1 consists of two semi - inf inite panels connected by a 
single edge at x - X s 0 (so that the summations of §2 reduce to one term, 
corresponding to N - 1, say). The acoustic pressure generated by the boundary 
layer is given by equation (2.8), where Green's function G is the sum (2.9) of 
the homogeneous plate Green’s function G q and the correction G s due to the 
panel edges. The contribution from G q is identical to the radiation produced 
by turbulent flow over a smooth plate, and will be ignored. In subsonic flows 
this is expected to be relatively weak compared to the sound generated by the 
interaction of the boundary layer pressures with the panel edges. The latter 
is determined by replacing G in (2.8) by G g . 


By making use of the representation (2.28) for a single edge we find 


P(x, t) 


where 


" (2?r) 2 p, 

c 0 l x l 


a 

?. I 


3(n,w) (u-Uk) — j -i w [t] 

S(k,r ,/e n ,a>)a, (u>) ( -k) J e L J 

T_ (k)JC_ (k,w) V " ' 0 3 j 


dkdo; 


(3.1) 


S(k k k u) - — i— f” S(y,r) e - 1{k l3' 1 +k !) , J- fk 3y3-"'l d s ydr, (3,2) 

1 2 (2ir) *J m<D 


is the space- time Fourier transform of the aeroacoustic source term, and [t] - 
t-|x|/c Q is the retarded time. The summation is over j - 0 to 3, where the 
shorthand notation a^ = A has been introduced; the minus signs in (-k) J and 
in the subscripts of , L_ have arisen because the integration variable k has 
been replaced by -k. 

The acoustic pressure can be related to empirical representations of the 
wall pressure fluctuations P w ( x 1 » x 3> t ) beneath the turbulent boundary layer by 
expressing S (k , T k n 3 ,co) in terms of the Fourier transform 

p w (k,/c 0 n 3 ,«) = (l/27r) 3 j^ i>i p w (x 1 ,x 3 1 t)exp{-i(kx 1 +/c 0 n 3 x 3 -wt)}dx 1 dx 3 dt 
of the pressure fluctuations that would be produced by the same turbulent 
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boundary layer flow over a smooth, rigid wall. The edge scattered pressure 
fluctuations (including sound and hydrodynamic pressures) then represent the 
correction of the wall pressure produced by the structural inhomogeneities. 


When the wall is rigid equations (2.2) and (2.3) are easily solved to 
supply 


Pw (k ’*0 n 3 


2nip 1 (w-Uk) 
a>r_ (k) 


S(k,T ,/c 0 n 3 ,w) . 


(3.3) 


This relation permits (3.1) to be cast in the form 


p(x, t) 


2ni 

c 0 l x l 


a 

S. J 


K^’ K o n 3 1 w)a.(w)(-k) J e- lw t t ^ dkdw 


1*1 


(3.4) 


which determines the interior noise pressure fluctuations in terms of the wall 
pressure fluctuations that are generated by the same turbulent flow on a rigid 
wall (the so-called boundary layer blocked pressure [10, 11, 18, 19]). 


The characteristic dimension of the boundary layer turbulence eddies is 
of the same order as the boundary layer displacement thickness 6* [20, 21]. 

When this is much smaller than the extent of the boundary layer the wall 
pressure p ( x ,x ,t) ma Y assumed to be statistically stationary, and to 
possess a wall pressure wavenumber- frequency spectrum ,k 3 ,u>) , which 

varies slowly with streamwise location x 1 . If a transverse section -|L < x 3 < 
|L of the elastic plate is wetted by the turbulent flow, where L > <5* , 

P (k, ,k„ ,«) satisfies 

r w 1 3 


(P W (k l lk 3 ’ W) Pw (K l ’ k 3 ^ <5(w-“’) <5 ( k 1 - K 1 ) P w ( k i. k 3’ W >> 6 * ^ Ll 


(3.5) 


The mean square acoustic pressure (p 2 (x,t)) in the interior may now be 

expressed in terms of a frequency spectrum density $(x,u>) (> 0) according to 

00 

(p 2 (x,t)) -* f $(x,u) do), (3.6) 

J o 

and equations (3.4) and (3.5), with account taken of the relations (2.29), 
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then yield 


$(x,w) 


4ttLw z |3(n,w) | 2 
C 2|x| 2 


a 

I 


P w (k,K 0 n 3 ,u>) 
\l_ (k,u>)|2 


llj a. (w)(-k) j | z dk, I x | 


(3.7) 


3.2 Clamped edge 

When the panel edges are clamped the displacement Z(y x ) of (2.26) must 
satisfy. Z = dZ/dy 1 = 0 as y ; -*• ±0. These conditions imply that Z and dZ/dy 1 
are continuous at the edge, and therefore that a 2 - a 3 s 0. Then a Q and a 1 
are determined by 

00 

a o I o + a i I i “ ' 1 ’ a o I i + a ! I 2 ” K o n i’ where - J k j dk/jf + (k,«) . (3.8) 


The solutions of these equations are conveniently expressed in the 
following dimensionless forms: 


nvco 2 



o 


mu> 2 

Q, , 

R2 1 


(3.9) 


_ *2 + *lM n l 

0 *2-^ >Ir 

1 0 2 


“l - 


-[* 1 + ♦o/<n 1 ] 

\T/2 _ \Tf \lr 
1 0 2 


where , 

A J dA 

*i = — 1 

._ w /, ie f 1 ^ (PoC§/p 1 c?)(c 0 /c 1+ MA/ / x)2 1 

U ' 1 *" -l - ^ [v(m 2 -A2) V{(c 0 /i/ Cl +MA)2-A 2 ) J 

A 2 =A 2 + M 2 n 2 , K = (mu 2 /D ) 1/4 > 0, M = V^T. (3.10) 

In these expressions K q is the in vacuo wavenumber of bending waves of radian 
frequency w. The value of the integrand is made unique by taking branch cuts 
for V(p 2 - A 2) from A - ±At(l-n 2 ) 1/2 to ±i® in the complex A-plane, and for 
V{ (c 0 ^/c 1 +MA) 2 -A 2 ) (when M < 1) from 

A - M [c 0 M/c 1 ±{(c 0 /c 1 ) 2 -(l-M 2 )n 2 } 1/2 ]/(l-M 2 ) to ±ioo; 
it may be assumed that u> > 0, and relations (2.29) can then be used to 
evaluate a Q and a 1 for w < 0. 
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3.3 Simply supported edge 

When the panel edges are simply supported the displacement satisfies Z 
d 2 Z/dy 2 = 0 as y 1 ■> ±0. In this case a 1 «= a 3 =0, and 


mo; 


mcj 


a ° K a o ' a 2 k 3 0=2 * 

o o 




a n — 

0 * ' 2 \Tf * 

1 0 2 


-[* 1 - »o(/ini) 2 ] 

\Tf * 2 _fl/> 1 

1 0 2 


(3.11) 


where , 


7TZ 


2j+l 


[i - Mr] 


- , 1 , (Picl/ppCg) (c 1 /c 0 +MA//i) 2 'j 2 j dX 

H lV (/i 2 -A 2 ) V{ (c 0 /i/c 1 +MA) 2 -A z ) J 


iif 

/H 




(3.12) 


and z = l/(l-irj) 1/4 . 

3.4 Numerical results 

To evaluate the integral in (3.7) for the acoustic pressure spectrum the 
functional form of the wall pressure wavenumber - frequency spectrum 
P w (k ,« 0 n 3 ,c °) must be specified. The dominant wall pressure fluctuations are 
associated with wavenumbers k in the convective domain, which peaks in the 
neighborhood of k - co/U c , where the mean convection velocity U c * 0.7U. We 
shall present numerical predictions of $(x,w) only for n 3 = x 3 /|x| - 0, i.e., 
for radiation directions in the plane of symmetry of the boundary layer flow. 
The crudest approximation for the dominant wall pressure fluctuations is then 
obtained by writing 


P w (k,0,u>) 



7T 


6(k- W /U c )* pp <«), 


(3.13) 
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where Z is a spanwise correlation length (in the x 3 -direction) of wall 
pressure fluctuations of frequency u>, and $ pp (to) is the point frequency 
spectrum, which satisfies <p^} - $ pp (w)dw. 


To 

case p Q 
$(x,w) 


further simplify the discussion, we shall limit consideration to the 
= p c o « c l . Then (3.7) becomes, 

(w/w c ) 3 |3(0,w) 


4Li„ M° 


|x|‘ 


I *IL 




eM 


V(l-M 2 ) 


eM J (M-M )“ 

c x c 

V[1-(M-M) 2 ] 


> 0, (3.14) 


where M = U /c is the convection Mach number, 0 = sin 1 (x /|x|) is the 

c c 9 0 1 

observer direction indicated in Figure 3, in the plane of symmetry x 3 ** 0, and 


3(0, w) 


ie 


[[d-ir,)^ ) sin 4 0-l]^ -ie 


1 ^ (1-Msin0) 2 

_cos0 + V{ (l-Msin0) 2 -sin 2 0) 


d) 


(3.15) 

Equation (3.14) expresses the acoustic pressure spectrum directly in 
terms of the turbulent wall pressure frequency spectrum $ pp (o>) . The ratio 


A(w) 


[*(x,to)/* pp («)] 

[ Lh/ | x | 2 ] 


(3.16) 


represents the efficiency of sound production by the interaction of the 
boundary layer with the panel edge. This is plotted as a function of in 

Figure 3 for a clamped edge connecting aluminum panels in air for three 
different mean flow Mach numbers M - U/c Q when 0 — 0, ty = 0.2 and the 
correlation length Z is given by the Corcos approximation Z 3 ** 1.4U c /w, U c — 
0.7U [10, 11, 18, 19, 22]. The efficiency decreases with increasing 
frequency, except at the resonance peaks, which occur where the convection 
velocity is close to the phase speed of flexural waves on the undamped plate; 
the magnitudes of these peaks are determined by the loss factor r } , whose 
influence is important only close to the resonance frequencies. The large 
value f} «= 0.2 used in these calculations is assumed to account for the 
relatively high damping of structural motions that would in practice be 
produced by the "trimming" of aircraft panels. The corresponding frequency 
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dependence of A for the simply supported edge is very similar to that shown in 
Figure 3 . 

To make quantitative predictions of the acoustic pressure spectrum $(x,o>) 
we use the following approximation for the boundary layer point pressure 
spectrum (w) proposed by Laganelli and Wolfe [12] (extending earlier 
analyses reported in [23, 24]) for an adiabatic wall: 

<“> . g2F ' 0 57 (3 17) 

(p o v 2)2 4 jt(v./U)M 1 + F 2 - 87 (w6*/U) 2 ] ’ 


where p Q is the mean density and v A is the friction velocity [20, 21] . The 
friction velocity varies with the boundary layer Reynolds number, but is 
typically of order 0.04U. Also, F = 1+0. 13M 2 , o = P Ims /\p\^ 2 * w here p rms is the 
root mean square wall pressure fluctuation. Recent studies indicate that o » 
0.01 [25], and this value is used below. Figures 4 and 5 depict predictions 
of the nondimens ional interior acoustic pressure spectrum 

*(x,w) = A(<j)$ pp (o>) 

O 0 2c *[Lh/|x| 2 ] 6*p*cl ’ 

respectively for clamped and simply supported edge conditions. The 
calculations are for a one mm thick aluminum plate at three different mean 
flow Mach numbers for 0 = 0°, r? ■= 0.2, when the boundary layer displacement 
thickness <5* is assumed to vary as 1/M 1/5 [20, 21], with 6* - 1 cm at M - 0.3. 
Both figures exhibit increasing acoustic power at low frequencies and 
resonance contributions at higher frequencies (the resonance being relatively 
more important at low Mach numbers) where the convection velocity U c coincides 
with the flexural wave speed. The overall levels are similar for either edge 
condition, especially at low frequencies where the coefficient a Q of (3.11) 
turns out to be very large compared to the other scattering coefficients a , 
j > 0. The predicted large sound power at low frequencies occurs because the 
Laganelli and Wolfe model spectrum (3.17) is effectively constant for 
o)6* /U < 1, which is consistent with an alternative empirical representation 
proposed by Efimtsov [26] (see also [27, 28]). Note, in particular, that the 
generated sound is significant only for o)6 * /U < 1, which is in accord with our 
assumption in §2. 
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4. PERIODIC PANELS WITH CLAMPED EDGES 

Consider next a plate divided into a periodic array of parallel, 
infinitely long panels of width d, with edges at x t = h Nd transverse to 
the mean flow direction (where N ranges over all integers in -® < N < <*>) . 

When the edge conditions are the same for all panels, periodicity implies that 
the coefficients A fIj in the representation (2.28) of the scattering Green’s 
function can be written 

A Nj = a je - iN \> n i d , j - 0, 1, 2, 3. (4.1) 

The equations determining the coefficients a. are obtained by application of 
the edge conditions at any one of the edges. When the edges are clamped Z «= 
<3Z/dy 1 = 0 at y 1 - Nd, and we then find that a 2 * a 3 =0, and that a Q , a 1 are 
given by (3.9) in terms of the following discrete analog of the definition 
(3.10) 


ief 1 (Pic?/> 0 co ) (co/ci+MAn/zx) 2 ) 

4V(m 2 -A2) V{(c 0 /x/ Ci +MA n )2-A2) J 

A n = 2N7r/K o d - /in, , A* = A 2 + /i 2 n 2 . (4.2) 

By proceeding as in §3.1 we may express the far field acoustic pressure 
in the interior region in the form 

( 2 tt) 2 i ” r" y to3(n,u>)(-l) J a J (o>) 

P X,t c 0 d l x l n-- L J l - (^ 0 n 1 +2Njr/d.«) 

x {« 0 n 1 +2N7r/d} J p w (n 1 +2N7r/d,/c 0 n 3 ,«)e’^ w ^ dw, |x| ■+ °°. (4.3) 

To use this result to calculate the acoustic frequency spectrum $(x,w) it 
is necessary to restrict the region of excitation of the panels by the 
boundary layer to a finite area A, since $ must increase linearly with A. 

When Va > d > 6* , there are contributions from many edges, and 


= 


2tt 

K~“d 


i 

N=-“ 
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(2*) 


2 HIT 


6 (w-w’ 



(4.4) 


where 6 is the Kronecker delta (equal to 1 for N = N’ , and 0 otherwise). 
The acoustic pressure frequency spectrum is then given by the following 
generalization of (3.7) 


8tt 2 Aw 2 |3(n,w) | 2 
c 2 d 2 | x | 2 


l 


N=~oo 


P M («Q n l + 2Hir/d,K 0 n 3 ,w) 

I £_ (^r^+ZNTr/d.w) | 2 




3j (w)(-l) j («; 0 n 1 +2N7r/d) j | 2 


(4.5) 

The value of the summation with respect to N depends on the functional 
form of the wall pressure wavenumber- frequency spectrum (k x , k 3 , w) . In the 

absence of precise measurements at Mach numbers of interest in the present 
application, we shall fall back on the simple model (3.13). This can be used 
in (4.5) provided it is assumed that the length scale of the turbulent motions 
is small compared to the distance d between the edges. The summation can then 
be approximated by an integral with respect to k, say, where we set k - 
/c n +2N?r/d # and dk - 27r/dL This is equivalent to assuming that the direct 
turbulence excitations at different edges are statistically independent. In 
the simplified case considered in §3, in which p Q - p , c Q - c , and n 3 = 0, 
we then obtain 


#(x,w) 


4Ai 3 M l(»/u> e ) 3 |3(0,«0| 2 |I j a J (^)(-M/M c ) J | 2 * pp (o>) 


d|x| 


>•*»>(§ T-'tfS • 


€M D 


eM* (M-M c ) 2 


C V(l-M‘) 


V[ 1 - (M-M ) ] 


a) > 0, (4.6) 


This is formally N p times larger than the corresponding expression (3.14) for 
a single, isolated panel edge, where N p - A/dL is the number of panel edges 
within the area A of the boundary layer. However, the functions (3.10) and 
(4.2) used in the definitions of the coefficients a Q and <* 1 are equivalent 
only when K Q d is large, i.e., when the characteristic bending wavelength of 
disturbances on a panel is small compared to the panel width, which occurs at 
high frequencies. At lower frequencies we might therefore expect predictions 
of (4.6) and (3.14) to differ substantially. 
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The efficiency A of sound generation can be defined by analogy with 
expression (3.16) for a single edge: 


A(w) 


[$(x,co)/$ pp (m)] 

[Ah/d|x| 2 ] 


(4.7) 


The variation of A(w) with w/co is depicted in Figure 6 for two mean flow Mach 
numbers M = 0.3 and 0.9 and 0 - 0° , for the aluminum plate considered 
previously in §3 (17 - 0.2) with panel edge separation of d - 50 cm. At very 
high frequencies, when flexural motions are rapidly damped, the efficiency of 
sound production at each edge would be expected to reduce to that shown in 
Figure 3 for an isolated edge. Interference between sound production by 
different edges causes the high frequency behavior to be oscillatory, and it 
is clear from the figures that, at these frequencies, it is only at the higher 
Mach number of M - 0.9 that the single and multiple edge radiation 
efficiencies are similar. The interference arises because flexural motions 
generated at a panel edge are subsequently scattered into sound at neighboring 
edges; alternatively, they may be interpreted in terms of resonance modes of 
the panels, since the peaks and troughs respectively occur where -flf is 
small and large, and where the coefficient a Q in (4.6) is large and small (o^/i 
is uniformly small). In all cases, however, the radiation efficiencies 
decrease very rapidly above the convective resonance frequencies indicated in 
the figure. At low frequencies the density of structural modes is small for 
multiple edges, and the efficiency is then very much smaller than for a single 
edge . 


The corresponding predictions of the interior acoustic pressure spectrum 

$(x,u) = A(w)$ pp («) 

6*pfc>[Ah/d|x|*] O 0 2c o * 

are shown in Figure 7 for M - 0.3, 0.6, 0.9, when the wall pressure frequency 
spectrum (w) is given by the Laganelli and Wolfe formula (3.17) (as in §3 
we assume the boundary layer displacement thickness 6 * varies like 1/M 1/5 , with 
6 * - 1 cm when M - 0.3). Comparison with the single edge case of Figure 4 
reveals that interference between edges causes large fluctuations in the 
spectrum about the single edge case, especially at lower Mach numbers; at low 
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frequencies the radiation is significantly smaller than for the single edge. 
The number of these fluctuations increases with increasing Mach number, their 
amplitudes being governed by the value of the loss factor rj , and are 
progressively larger as rj -* 0. At very high frequencies the amplitude of the 
sound, although fluctuating rapidly, does not vary significantly with Mach 
number. As before, it should be noted that significant generation of sound is 
predicted only for /U < 1. 

The directivity of the sound is illustrated in Figure 8 for w<5*/U - 0.05. 
These plots have been made by using the Laganelli and Wolfe spectrum (3.17) at 
a constant displacement thickness 6* = 1 cm. The directivity exhibits a 
uniform, monopole type of behavior in the plane of symmetry of the flow (x 3 - 
0) , except for a very narrow null at 0 = sin 1 (1/(1+M)). This angle 
corresponds to the angle of incidence in the reciprocal, reverse flow problem 
at which the pressure differential [p R ], say, across the homogeneous plate 
vanishes identically, and the associated flexural response $" 0 of equation 
(2.16) is null. These results are, of course, strictly relevant only in the 
flat plate approximation. 
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5. CONCLUSION 

An important component of mid- to -high frequency aircraft interior noise 
is produced by the scattering of intense boundary layer wall pressures by 
structural inhomogeneities of the fuselage. In this chapter an idealized 
configuration has been investigated in which the fuselage is plane, and formed 
by thin elastic panels of infinite length in the direction transverse to the 
mean flow. This may be regarded as a local model describing the influence of 
fuselage ring stiffeners, which usually constitute the most significant 
constraint on panel motions. The excitation of the panels by the high 
Reynolds number boundary layer flow was expressed in terms of a prescribed 
boundary layer wall pressure spectrum. The sound produced by the interaction 
of these wall pressures with the panels and their edges has been determined 
analytically for the case where all panels have equal width and are subject to 
the same edge conditions. The procedure is expected to be valid provided the 
principal acoustic wavelengths are larger than the boundary layer displacment 
thickness, which is the case at subsonic mean velocities provided wd /U < 1 , 
and the detailed predictions indicate that sound generation at higher 
frequencies is negligible. The numerical results give explicit 
representations of the influence of the main stream Mach number on sound 
generation. Detailed results of this type can be used to validate more 
general numerical prediction schemes . 

The material of this chpater is an expanded version of a presentation 
made at an Interior Noise Workshop organized by Dr. R. J. Silcox at NASA 
Langley Research Center (25 - 27 April, 1995). 
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APPENDIX: REVERSE FLOW RECIPROCITY 


For given points x A and x^ respectively in the interior (x 2 > 0) and 
exterior (x 2 < 0) domains of Figure 2, the reciprocal theorem (2.13) can be 
stated in the form 

G(x a ,x b ,w) - G R (x fi ,x A ,w) , (A.l) 

p o 

where G and G R represent perturbation stagnation enthalpies associated 
respectively with the direct problem 


[V 2 + « 2 ] G (x , x B ,w) =0, 

[V 2 + (k 1 + iMa/dx 1 ) 2 JGCx.Xg ,u) = dCx-Xg) 


x 2 > 0, 


X 2 • X B2 < 0 ■ 


(A. 2) 


and the reciprocal (or adjoint) problem 


[V 2 + /c 2 ]G R (x,x w) = 6 (x-x ) , x , x > 0, 


[V 2 + (k 1 - iMa/5x 1 y ]G r (x,x a ,w) “ 0, 


x 2 , <0. 


(A. 3) 


Let p denote the local mean fluid density. In the usual way [29], 
multiply equations (A. 2) by pG R (x,x A ,u>) and (A. 3) by pG(x,x fi ,cj) , subtract 
corresponding resulting equations in (A. 2) and (A. 3), and integrate with 
respect to x over their respective regions of validity. The divergence 
theorem is then applied to transform the volume integrals into integrals over 
the upper and lower surfaces of the plate (the radiation condition ensures 
that there are no contributions from surfaces at infinity). 


If f and denote the flexural displacements of the plate in the direct 
and reciprocal problems, then 

8G/dx 2 - -w 2 f, dG R /ax 2 - -w 2 r R , x 2 - +0, 

dG/dx 2 - iw( -iw+Ua/3x 1 )f , 9G R /3x 2 - iw( - iw-U3/3x 1 ) f R , x 2 •* -0. 

These relations are used to eliminate normal derivatives of G and G R from the 
surface integral. By integrating by parts with respect to x x and making use 
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of the time -harmonic form of (2.3) we then find 

-w 2 J [r [p R ] - r R [p])dx 1 dx 3 - P 1 G r (x b ,x a ,w) - P 0 G(x a ,x b ,w), (A. 4) 

- 00 

where the integration is over the undisturbed plane x £ - 0 of the plate, and 
[p], [P R ] are defined as in (2.5). Equation (A.l) is now deduced by 
substituting for [p] and [p R ] from the corresponding time-harmonic bending 
wave equation (2.4), and observing that the resulting integral vanishes 
identically. In doing this, note that contributions to the integral in (A. 4) 
from the immediate neighborhoods of the edges are null since the displacements 
and pressure discontinuities are finite. 
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Figure 1. Configuration of the turbulent flow and inhomogeneous elastic 
plate . 
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Figure 2 . 
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Figure 4. Predicted interior sound pressure spectrum level 

10xLg[4>(x, w )/{6*p2c3Lh/|x|2}] 

for an isolated clamped aluminum panel edge in air when 0 = 0° , 
h - 0.1 cm, tf ^ 0.2. The displacement thickness 6* is assumed 
to vary as 1/M 1/5 , with 6* - 1 cm when M -* 0.3. 
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Figure 5. Predicted interior sound pressure spectrum level 

10xLg[$(x, W )/{6*p*c*Lh/|x|2}] 

for an isolated simply supported aluminum panel edge in air 
when 0 = 0°, h - 0.1 cm, rf - 0.2. The displacement thickness 
6* is assumed to vary as 1/M 1/5 , with 6* » 1 cm when M - 0.3. 
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Figure 7. Predicted interior sound pressure spectrum level 

10xLg[$(x, w)/{ 6* pjj Ah/d|x| 2 } ] 

for a periodically clamped aluminum plate in air when 0-0°, 
h = 0.1 cm, d = 50 cm, »? - 0.2, and when $> pp (w) is given by the 
Laganelli and Wolfe formula (3.17). The displacement thickness 
6* is assumed to vary as 1/M 1/5 , with 6* - 1 cm when M - 0.3. 
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Figure 8. Predicted interior noise directivity (for x 3 - 0) for a 

periodically clamped aluminum plate in air when h - 0.1 cm, 
d - 50 cm, r) — 0.2, 6* - 1 cm, and * pp (w) is given by (3.17). 


36 




Report No. AM-97-027 


Boston University, College of Engineering 


CHAPTER 2 

SOUND GENERATED BY A VORTEX INTERACTING 
WITH A RIB-STIFFENED ELASTIC PLATE 
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SUMMARY 

An analytical investigation is made of the sound produced when a line 
vortex translates over a parallel, nominally smooth rib-stiffener on a thin 
elastic plate. This is a canonical fluid-structure interaction that is 
believed to be an important source of boundary layer generated aircraft 
interior noise. To eliminate additional complications introduced by mean 
flow, the vortex motion is assumed to be controlled by its image in the plate. 
The arrival of the vortex at the rib is preceded by bending wave forerunners 
which are excited in the plate at that frequency Q at which the flexural wave 
phase velocity is just equal to the vortex translational velocity v. These 
waves have group velocity equal to 2u, and extend over a length of plate ahead 
of the vortex determined by the structural damping. The principal acoustic 
source is the interaction of this forerunner with the stiffener, and the sound 
is predominantly of frequency Q. The amplitude of the radiation progressively 
increases to a maximum as the vortex approaches the rib, and subsequently 
decreases rapidly to zero with passage of the vortex over the rib. Numerical 
results for the efficiency of sound generation, its directivity, and the 
acoustic pressure signature, are given for clamped and simply supported 
conditions at the stiffener for aluminum plates in air. 
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1 . INTRODUCTION 

A variety of noise sources contribute to aircraft cabin noise , and are 
generally classified as airborne or structure -borne sources [1]. Structure- 
borne sound is caused by engine vibration and by vibrations induced by 
impingement of engine exhausts on wings and flow control surfaces, and is 
transmitted to the fuselage via vibrational motions of the airframe. Airborne 
sources -include the direct incidence of engine generated sound on cabin walls 
and the high speed turbulent air flow over the exterior fuselage. A critical 
f olg is played by the frequency-dependent transmission characteristics of the 
fuselage which are manifested by irregular variations of sound levels with 
frequency and location within the cabin [2]. Fuselage ring- stiffeners , 
stringers, windows, etc, have a significant influence on the modal response of 
the cabin to outside forcing and incident sound, but they can also behave as 
sources of sound when they interact with the exterior turbulent flow and with 
structural vibrations induced by the flow. This boundary layer generated 
noise is believed to make a substantial contribution to the interior noise of 
a jet transport in steady cruise, and is currently the subject of intense 
experimental scrutiny and numerical modeling [3 - 7]. 

Numerical schemes for predicting boundary layer generated interior noise 
must take proper account not only of sound waves produced by the fluid- 
structure interaction, but also of structural vibrations generated by the flow 
that can subsequently generate sound by their interactions with structural 
inhomogeneities. Such predictions must ultimately be validated by comparison 
with experiment, but it is also useful to have available a collection of 
benchmark analytic solutions against which the numerical results can be 
compared . 

In this chapter we investigate a canonical problem that includes all of 
the major aspects of the fluid-structure interactions involved in boundary 
layer generated interior noise, yet is simple enough to be treated 
analytically. This problem involves an idealized boundary layer disturbance 
in the form of a line vortex. The vortex is adjacent to a section of the 
fuselage which is taken to be a locally flat elastic plate, and is propelled 
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under the influence of its image in the plate towards a rib stiffener where 
the plate is either clamped or simply supported. The stiffener is assumed to 
be parallel to the vortex, which further simplifies the analysis of their 
interaction, which is then governed by two-dimensional equations of motion. 

In a first approximation the vortex translates at a constant speed u, say, 
over the plate and generates predominantly bending waves of a fixed frequency 
whose phase velocity is u. However, for an elastic plate in air the group 
velocity of these waves is about 2u, so that the arrival of the vortex at the 
stiffener is preceded by flexural vibrations of fixed frequency and wavelength 
which generate sound by their interaction with the stiffener. We calculate 
the amplitude and waveform of the scattered sound, and determine the 
dependence of the efficiency of sound generation and the acoustic field shape 
on frequency. 

The fluid-structure interaction problem is formulated and solved for an 
arbitrary two-dimensional flow inhomogeneity adjacent to a rib- stiffened 
infinite plate in §2. In §3 we discuss the frequency dependence of the 
radiated sound for clamped and simply supported conditions at the stiffener , 
and consider in §4 specific results related to a vortex convecting at subsonic 
velocity over the stiffener. 
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2. EQUATIONS OF MOTION 


2.1 Formulation of the problem 


Consider the generation of sound by aeroacoustic sources adjacent to a 
thin, elastic plate which coincides with the plane x £ ■= 0 of the rectangular 
coordinate system (x^ ,x 2 ,x^ ) . The fluid on both sides of the plate is at rest 
in the undisturbed state, with mean density and sound speed respectively equal 
to p Q and c Q (Figure 1). The plate is divided into two semi - inf inite , 
homogeneous "panels" which are smoothly joined by a rib coinciding with the 
x^-axis, where their abutting edges are assumed to be either clamped or simply 
supported. It is required to determine the influence of this joint on the 
production of sound. 

To fix ideas, suppose the aeroacoustic sources are confined to the region 
x 2 > 0, and, with t denoting time, let p w (x 1 ,x 3 ,t) be the blocked pressure [8 
- 12], i.e., the pressure that these sources would exert on the surface x 2 - 
+0 of the plate when the latter is assumed to be rigid. In principle this 
pressure distribution can always be calculated, both on the plate and within 
the fluid in x £ >0, when the sources are known, which we shall assume to be 
the case. Since the blocked pressure is the solution of the governing 
equations in the absence of surface motion, it may be taken as a prescribed 
surface driving stress in the actual problem involving the ribbed, elastic 
plate . 


When the surface motion is only linearly disturbed from its planar form 
by the blocked pressure, the sound produced by the fluid-structure interaction 
can be determined at an arbitrary point x in the fluid by application of 
Rayleigh’s reciprocal theorem [13, 14]. Write 

P« (x i- x 3’ t) “ J*-Pw( x i ' x 3 • w > e ’ U ' t du> < (1) 

and consider a reciprocal problem in which a time harmonic volume point source 

q R e -iut is pq ace d at x. Then Rayleigh’s theorem asserts that 


P (x, w)q R - -iwj p w (y x ,y 3 ,y 3 ,w)dy 1 dy 3 , 


( 2 ) 


where f R (y, ,y 3 , w)e" iut is the displacement of the plate (in the x 2 -direction) 
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at the point (y ,y 3 ) produced by the source q R , and p(x,«)e is the acoustic 
pressure generated by the harmonic component P w ( x 1 > x 3 > w ) e iut of the blocked 
pressure . 


2.2 The reciprocal problem in two dimensions 


We now specialize the discussion to two-dimensional source distributions, 
which are independent of the coordinate x 3 parallel to the rib, so that here 
and henceforth x - (x 1 ,x 2 ), etc. The principal source of sound attributable 
to the fluid-structure interaction is confined to a neighborhood of the rib, 
near the coordinate origin, and equation (2) will be applied to determine the 
radiation at large distances. The propagation of sound of frequency w is 
governed by the Helmholtz equation 

(V 2 + k 2 )p - 0, K 0 - w/c 0 , (3) 


where the acoustic wavenumber k q may be assumed to have a small positive 
imaginary part. Small amplitude, time harmonic motion of the plate satisfies 
the thin plate bending wave equation [14, 15] 


(BaV^xJ + mw 2 )f + [p] =0, x 1 * 0, (4) 

where B is the bending stiffness, m the mass per unit area of the plate, f the 
flexural displacement of the plate (in the x 2 direction), and 

[p] - p^.+O.w) - p(x x , -0 ,w) (5) 

is the pressure loading. To account for damping in the plate we write 

B - B {1 - ir?(w) ) , (6) 

O 

where B is real and r?(w) is a loss factor satisfying tj(-u) - -rj(to) which will 

O 

be assumed to be of constant absolute value. 


In the reciprocal problem, q R is taken to be a line source of unit 
strength located at x. The velocity potential produced by this source at the 
point y in an unbounded medium is 

<pj(y,v) - -X^olx-yl). 

where is a Hankel function [16]. As the observer position x recedes 

to infinity, and for y near the plate, 
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- <P 0 e' i,c o n *y , 


. p i(* 0 |x|+ir/4) 


<P n 


V(8tt«; | x| ) 


1*1 


( 8 ) 


where n = (sin0,cos0) and 0 is the observer angle shown in Figure 1. We shall 
assume that x 2 > 0; the case of an observer in the lower region x 2 < 0 is 
obtained by changing the sign of n 2 . Scattered waves are produced when the 
plane wave (8) impinges on the plate. When the influence of the rib is 
ignored/ the scattered waves are plane reflected and transmitted waves, and 
the reciprocal velocity potential near the plate would then be 


<P 


R 

o 


<P o [e- iK o n 2 y 2 + Re i/c o n 2 ^ 2 ] e " i/c o n i , y 2 > 0, 

Te- i *o< n iyi +n 2 y 2 > , y. < o, 

o * 


(9) 


where R and T are suitable reflection and transmission coefficients. These 
coefficients are found in the usual way [15], by making use of the bending 
wave equation (4), and the linear relations p *= ip Q oxp t -icof = dcp/d y 2 , between 
pressure, velocity potential and displacement at y 2 = ±0, leading to 


T(«,0> 


1 - R(u>,0) 


-2ie 

{ (1-ir;) (w/o) c ) 2 sin^0 - 1 } (w/ci> c ) cos0 - 2ie’ 


(10) 


where u> c = c 2 (m/B Q ) 1/2 is the coincidence frequency of the undamped plate [15], 
above which the phase velocity of flexural waves on the plate in vacuo exceeds 
the speed of sound c Q , and e = P 0 C 0 A> C m is a fluid loading parameter that 
depends only on the material properties of the fluid and plate [17]. By 
making use of the general relation - iwf - d(p/d y 2 (y 2 ■* ±0), we find that the 
reciprocal displacement of the plate when the rib is ignored is given by 


r: - 


ri e- iK o n iyi 


-n 2 T(a,,0)e i < ,c o l x l +,r / 4) 

c 0 V(8?r/c 0 |x|) 


(ID 


Because of the rib, the net reciprocal displacement , say, of the plate 
includes an additional contribution that may be interpreted as the flexural 
motions induced by the scattering of f* at the rib. The general solution of 
scattering problems of this type is discussed in [18, 19], where it is shown 
that 
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cVr, 


oo 

I — 


k“e lk yi dk 
D(k,w) 


+ e 


-i/c 


o n iyi 


( 12 ) 


where the summation is from n - 0 to 3. The four "scattering" coefficients a n 
are determined by the mechanical conditions at the rib. They are independent 
of k and satisfy a n (-w) - (-l) n a*(w), where the asterisk denotes complex 
conjugate (see [6]). D(k,w) is the dispersion function for flexural waves on 
a fluid loaded plate [15]: 

D(k,w) s Bk* - nw 2 - 2 P ^J-/{k\ -k 2 ) 1/z . (13) 

The sound received at x as a result of the fluid- structure interaction is 
now obtained by substituting from (12) into Rayleigh’s formula (2). The total 
radiation received at x is found by adding to this prediction the radiation 
associated with the interaction of the aeroacoustic sources with a rigid wall, 
i.e, by adding the acoustic far field of the blocked pressure p w . Thus, the 
resulting expression for the acoustic pressure consists of two distinct terms: 
the first involves p w and the component of (2) due to the flexural 
displacement f* (corresponding to the second term on the right of equation 
(12)). This is just the sound produced by the same aeroacoustic sources 
adjacent to a smooth, homogeneous plate with no rib. For this component of 
the radiation the plate behaves as a passive reflector of the sound that would 
be generated by the same sources if the plate were removed. This is because 
there is no "wavenumber conversion" mechanism at a plane boundary, i.e., the 
energy of an incident disturbance of given trace wavenumber remains in that 
wavenumber after interacting with the plate, and only incident energy that is 
already sound will be re-radiated from the plate as sound. In particular, 
this radiation is null when the source consists of a steady disturbance 
translating over the plate at a constant subsonic velocity. Our applications 
in §§3,4 will be to sources of this type, and we shall therefore discard this 
term. 


The remaining contribution to the scattered sound is determined by the 
reciprocal flexural displacement generated at the rib, i.e., by the first term 
on the right of (12), and represents the interaction of the aeroacoustic 
sources with the rib. By substituting from (12) into (2), retaining only this 
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contribution, and integrating over all frequency w, the fluid- structure 
interaction noise is accordingly found to be given by 


oo 

7rin 5 f ~ 

P(X,t) = V(2 tt|x|) J ffi 


Vi/c^"T (w , 0)a n ( -k) n p w (k, w)e* iw t ^ 


D(k,w) 


dkdw , | x | -* ® , 


(14) 


where [t] = t-|x|/c is the retarded time, and 


00 

p w (k,co) - (1/2 tt) 2 f p w ( Xl ,t)e' l(kx i" ut) dx^t, 
-00 


is the space-time Fourier transform of the blocked pressure. 


(15) 


2.3 Scattering coefficients for clamped and simply supported rib conditions 


The coefficients a are determined from the mechanical conditions at the 

n 

rib. For clamped and simply supported conditions f R satisfies respectively 
- df R /3y 1 = 0 and f R - d 2 ^/dy\ = 0 as y x ±0 . 

Introduce dimensionless coefficients a defined by 

n J 

a = ( mw 2 /K n+1 ) a , (16) 

where K q = (mw 2 /!^) 1 ^ > 0 is the in vacuo wavenumber of bending waves of 
radian frequency u>. The evaluation of the a using the representation (12) 
can be accomplished by the method of reference [18], which supplies: 


Clamped 

“o = -VV “l = '"VV tt 2 “ °3 ■ °* 

where , 


n 


Jo ( 1 - ir; ) A 4 


A n dA 

1 -2i€//iV(/i 2 - A 2 ) * 




(17) 


(18) 


It may be assumed that a? > 0 in applying these formulae; the values of a^ = 
a^ (w) for w < 0 are obtained from the relation a n (-u>) « (-l) n a*(w). The 
integration contour in (18) runs below all singularities of the integrand on 
the real axis; in particular, the branch of V(/i 2 -A 2 ) which is positive 
imaginary for A > p (> 0) is taken. 
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Simply supported 


*■- 

2 

1 0 2 


Q. - 


[*i- 1 


tfr * 2 . 
1 


\Jr> 
0 2 




- 0 , 


where , 

wz 2n+1 n r 00 a 2 " dx 

= r i - C-l) 1 + 4ie — 

2 1 V ' J Jo { [ Cl-ir?)A 4 -1 ] M V(m 2 -A2) - 2ie}{ ( 1 - ir/) A 4 -1} 

and z - l/(l-ir ? ) 1/ \ 


(19) 


( 20 ) 
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3. SOUND GENERATED BY A TIME HARMONIC GUST 


The blocked pressure generated by an arbitrary, uniformly convecting 
aerodynamic source, which is independent of x 3 , is equivalent to a 
superposition of elementary time harmonic "gusts" of the form 


p ( W )e' iw < t - X i/ U) 


( 21 ) 


where v - constant is the convection velocity, which is assumed to be subsonic 
(< c Q ) . In §4 we shall use such a decomposition for a source consisting of a 
line vortex. Here we shall investigate the sound produced by the individual 
Fourier components (21) and their dependence on u. 


It follows immediately from equation (14) that the far field sound 
produced by the gust (21) is given by 


P(x, t) 

p o V(h/|x| ) 


- V( 7t/ 2 ) cosQT(w, 0)M i| p 2 (-p/ M) n c* n e + 

V(w c h/c 0 ){[(l-ir,)/-M*]M 2 - 2eM 5 //(l-M z ) } 


M = Vw/w. , 

( 22 ) 


where h is the thickness of the plate, and M = i>/c Q is the convection Mach 
number . 


The dependence of p(x,t) on the radiation direction 0 is governed by the 
transmission coefficient T(u;,0) (defined in (10)) and the ©-dependence of the 
scattering coefficients . The influence of the latter turns out to be 
small, and equations (10) and (22) therefore imply that the radiation will 
tend to have a dipole character, with p a cos0, when fx 2 s <V W C < € » whereas at 
higher frequencies (but less than cj c ) the field shape resembles that of a 
monopole source. These conclusions are illustrated in Figure 2 for an 
aluminum plate in air, for which 

u, c h/ c o - 0.22, e - 0.0021. (23) 

Thus, dipole behavior is expected for v>/v> < 10~ 3 . In Figure 2 the plate is 
assumed to be clamped at the rib, and |p(x,t)| z is plotted as a function of 0 
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for a convection Mach number M - 0.5, although the actual dependence on Mach 
number is small. The plots are scaled to the same level in the normal direc- 
tion 0 * 0°. The directivities for simply supported conditions at the rib 
differ imperceptibly from those shown in Figure 2. At very high frequencies 
(not shown) , when u> > o> c , the directivity is ultimately dominated by the 
radiation from "leaky" flexural waves generated at the rib. These waves have 
supersonic phase velocities and decay by radiating their energy into the fluid 
as sound. They correspond to complex roots of the denominator of expression 
(10) for T(u>,0) that occur near the real angles 0 - tsin" 1 {V(u> c /u>) } when 
< 1 , so that the radiation pattern is consists of two lobes centered on 
0 * isin' 1 [V(u) c /u) } . However, the efficiency with which sound is generated at 
these frequencies is very small, since when w > u> c flexural motions of the 
plate induced by the gust become negligible. 

The efficiency of sound production at different frequencies is 
illustrated in Figures 3 and 4, where 10xlog{ | p(x, t)/p o | 2 /(h/ | x | ) } (dB) is 
plotted against for 0=0°, and for three different convection Mach 

numbers, respectively for clamped and simply supported conditions. The 
"resonance" peaks occur near u>/w c * M 2 , where the real part of the denominator 
on the right of (22) vanishes; at this frequency the gust convection velocity 
v coincides with the bending wave phase speed. The magnitudes of the peaks 
depend on the loss factor rj which has been assigned a nominal value of 0.01. 

At higher frequencies the efficiency of sound production by scattering 
decreases precipitously. 
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4 . VORTEX SOUND 


Turn attention now to the special case in which the gust of Figure 1 is a 
line vortex of circulation I\ The vortex is assumed to translate over the 
plate with velocity induced by images in the plate. If the plate were rigid 
the translational velocity would be v — r/47rd, where d is the distance of the 
vortex from the plate, which would then be constant. The order of magnitude 
of the displacement of the plate beneath the vortex can be estimated from the 
bending wave equation (4) to be given by 

f/ d » (d/h) p Q xj z /[p s c z (h/d) 2 + p s v z ] (24) 

where p s is the density of the plate material, and c is the velocity of 
longitudinal plate waves. In aircraft applications the distance d may be 
identified with a boundary layer displacement thickness, and the ratio d/h is 
then typically of order 10. For an aluminum plate in air and a vortex 
velocity u < c Q , equation (24) then implies that f/d < 1. In a first 
approximation it is therefore reasonable to neglect the compliance of the 
plate in calculating the self- induced motion of the vortex, at least during 
the interval of time in which the vortex interacts significantly with the rib. 


In these circumstances the vortex may be imagined to exert a steady 
blocked pressure distribution on the plate that translates uniformly at speed 
v = r/47rd. We shall assume that the vortex translates in the positive 
x 1 - direction , and passes over the rib at time t - 0. The blocked pressure 
field of the vortex will be approximated by that for a vortex in 
incompressible flow; any differences due to the effects of finite Mach number 
will modify the overall level of the acoustic radiation, but should not affect 
its general characteristics. Thus [20] 


P w ( x i.t) - p o” 2 


4d 2 {(x 1 -ut) 2 -d 2 } 
{ (x x -ut) 2 +d 2 } 2 


and 


p w (k,w) = -2p 0 d 2 | w| e’ l ud / u 1 6(k-w/u) . 


(25) 
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Before proceeding to the evaluation of the radiation integral (14), 
further insight into the expected form of the interaction radiation can be 
obtained by considering the flexural wave motions produced by vortex. These 
waves generate sound by their subsequent impingement on the rib. To do this, 
replace [p] by [p] + p w in the bending wave equation (4), where p is a 
solution of equation (3) with "outgoing" wave behavior. When the presence of 
the rib is ignored, the flexural displacement produced by the vortex is 
then found to be 


f 


W 



|fa)|e -|«d/u|-io(t-x 1 /w) ^ 
D(w/u , w) 


( 26 ) 


At large distances from the vortex the value of this integral is dominated by 
residue contributions from poles of the integrand near the real w-axis. These 
poles are located above the real axis at w = ±u 2 {m/[B o (lTirj) ] } 1/2 , provided the 
small influence on their locations of fluid loading is neglected. It follows 
that 

f w - f° sin{£2(t-x x /v) }exp{ -|»jQ(x 1 /u- 1) ) , x 1 -ut -* +® 

= 0, x. -ut ■* -°°, (27) 

where £2 = v z (m/B ) 1/2 = M 2 w c (M = v/c Q ) , and is a constant determined by the 
residues at the poles. 

The frequency £2 is that of the flexural wave resonantly excited by the 
vortex, whose phase velocity is the same as the vortex translation velocity u. 
However, the group velocity [21, 22] 

du/dk = -OD(k,w)/0k)/(3D(k,w)/3w) » 2k(B o /m) 1/2 
of this wave is 2u, so that energy imparted to the surface by the vortex runs 
ahead of the vortex. In the absence of damping (»j ■* 0) this flexural wave 
"forerunner" would extend infinitely far ahead of the current location Xj^ - ut 
of the vortex, and its interaction with the rib would generate time harmonic 
radiation of constant strength and frequency £2 for all times prior to the 
arrival of the vortex at the rib. According to (27), the flexural wave energy 
in the wake of the vortex is negligible, so that the acoustic radiation should 
decay rapidly after the passage of the vortex over the rib. It should be 
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noted that this conclusion is in no way an artifact of the present idealized 
model; energy supplied to a plate by boundary layer excitation at hydrodynamic 
coincidence [11] will also tend to propagate ahead of the convecting forcing 
region, although this is not usually evident when the interaction is analyzed 
solely in the wavenumber- frequency domain. 

These qualitative predictions are borne out by numerical evaluation of 
the radiation integral (14), with p w (k,«) given by (25). The integral can be 
expressed in the following form 


p(x, t) 


P 0 c £V(h/|x| ) 


Re^/(0,M,z)e" lzu t / d j dz, 


where M = v/c , 


/(0,M,z) = 


2V(27r)M 6 cos0T(w,0)zp 2 X n (*M/M) n a n e' (z+1,r/4) 


V(a> c h/c 0 )[(l-iT7)M b 


M 4 m 2 


2eM 5 /(l-M 2 ) 1/2 ] ’ 


(28) 


(29) 


and p - Vc JfiT, w/w c - zM(h/d)/(w c h/c Q ) . Observe that /(0,M,z) has a pole 
close to the real axis near w/<*> c - M 2 when e < 1, i.e. , at the resonance 
frequency fi. 

Figures 5 and 6 depict calculated acoustic pressure signatures (28) for 
clamped and simply supported rib conditions at several subsonic Mach numbers M 
- u/c Q . The calculations have been performed for an aluminum plate in air for 

h/d - 0.1, r] - 0.05 and 0-0°, 

and the curves are plotted as a function of the retarded position of the 
vortex scaled on the stand-off distance d. The acoustic pressure is dominated 
by the resonant mode of frequency Q excited by the vortex. In terms of the 
efficiency plots in Figures 3 and 4, this means that the dominant radiation is 
from the neighborhood of the peak at w/u> c - M 2 . The acoustic amplitude grows 
exponentially with time until it is cut-off when the vortex passes over the 
rib at t - 0. The growth rate of the sound is determined by the loss factor 
Tj] a value of 0.05 is large for a homogeneous aluminum plate, but is probably 
representative of damped aircraft structures, and has been used to illustrate 
the role of the flexural wave forerunner (27) ahead of the vortex. When rj is 
fixed the damping of the forerunner increases with frequency, so that the low 
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Mach number (M - 0.1) radiation shown in the figures grows much more slowly 
than the higher Mach number waveforms. The characteristic acoustic wavelength 
* 27rc Q /fi = 27rd/[M z (w c h/c Q ) (d/h) ] - 285d, 32d, 8d respectively for the three 
Mach numbers M - 0.1, 0.3 and 0.6 shown in the figures, i.e., the wavelength 
of the dominant radiation is typically large compared to the stand-off 
distance d of the vortex. The plots in Figures 5 and 6 are drawn to the same 
nominal scale. The maximum acoustic amplitude therefore increases with Mach 
number, and the amplitude of the sound is marginally smaller for the simply 
supported rib condition, which is in accord with the efficiency plots of 
Figures 3 and 4. The directional characteristics of the sound can be surmised 
from Figure 2 by noting that £2/w c “ 10' 2 , 10' 1 , 0.36 respectively for M - 0.1, 
0.3 and 0.6. Thus, higher Mach number radiation will exhibit a monopole field 
shape; at very low Mach numbers a dipole directivity would be expected. 
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5. CONCLUSION 

The sound produced by turbulent flow over the fuselage of an aircraft is 
a significant component of aircraft cabin interior noise. In this chapter we 
have examined a simple canonical problem that models an important aspect of 
this fluid-structure interaction, in which a hydrodynamic disturbance in the 
form of a line vortex convects over a parallel rib on an elastic plate. 

The arrival of the vortex at the rib is preceded by a flexural wave 
forerunner in the plate of frequency Q « M 2 u> c , whose phase velocity is equal 
to the vortex convection velocity, M being the convection Mach number and w c 
the coincidence frequency. The group velocity of this structural wave is 
twice the vortex translation velocity, and the forerunner extends a distance 
ahead of the vortex determined by the structural damping. There is no 
preferred frequency associated with the vortex motion, but the interaction 
radiation is dominated by sound of frequency Q, whose amplitude increases 
exponentially as the vortex approaches the rib, but is subsequently rapidly 
annulled after passing the rib. Detailed predictions have been worked out for 
clamped and simply supported rib conditions for an aluminum plate in air, and 
in both cases indicate that the characteristics of the sound are dominated by 
the condition that the plate displacement is required to vanish at the rib. 
Conditions on the slope and curvature of the plate at the rib do not 
appreciably affect the radiation, although the radiation from a clamped rib is 
larger by a few dB. The directivity of the sound depends on the value of 
Q/ct> c , i.e. on the translational Mach number of the vortex; it is effectively 
the same for clamped and simply supported joints. 
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x 2 



Figure 1. Aeroacoustic sources adjacent to a rib-stiffened thin elastic 
plate . 
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rib 

Aluminum/ Air 


Figure 2. Directivity of the scattered sound for clamped rib conditions: 
aluminum/air, for a convection Mach number M - 0.5 and 
w/w c = 10' 1 , 10' 2 , 10* 3 . 
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Figure 3. Efficiency of sound generation for clamped rib conditions: 
aluminum/air, for tj - 0.01 and M - 0.2, 0.5, 0.8. 
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Figure 4. Efficiency of sound generation for simply supported rib 
conditions for 0 = 0°: aluminum/air, for rj - 0.01 
and M - 0.2, 0.5, 0.8. 
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Figure 5. Acoustic pressure signatures (to the same nominal scale) of sound 
generated for clamped rib conditions, plotted as a function of the 
nondimensional retarded position u[t]/d of the vortex for different 
translation Mach numbers M - u/c Q , and for 0 « 0°, h/d *0.1 
and r) = 0.05. The vortex passes over the rib at t = 0 . 
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Figure 6. Acoustic pressure signatures (to the same nominal scale) of sound 

generated for simply supported conditions, plotted as a function of 
the nondimensional retarded position u[t]/d of the vortex for 
different translation Mach numbers M — u/c Q , and for 0-0°, 
h/d = 0.1 and rj = 0.05. 
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CHAPTER 3 

ON THE CONTRIBUTION FROM SKIN STEPS TO 
BOUNDARY LAYER GENERATED INTERIOR NOISE 
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SUMMARY 

An estimate is made of the aircraft interior noise produced by high subsonic turbulent 
wall pressures interacting with a fuselage skin step formed when adjacent, elastic panels 
overlap. The panels are modeled as thin elastic plates bonded in the vicinity of the step, and 
clamped along a line transverse to the mean flow direction. Sound is produced by scattering 
of the turbulence induced flexural skin motions at the clamp line, and by interaction of 
the turbulence pressure with the step. The skin step component is found to be significant 
at high frequencies, above the convective resonance frequency of the panels (at which the 
turbulence convection velocity equals the flexural wave speed). At typical cruise Mach 
numbers, however, when boundary layer generated sound is believed to dominate interior 
cabin noise, the critical frequency above which skin step noise is important is of order 10 
kHz, which is beyond the range of interest in practice. The critical frequency is much lower 
at take-off and landing approach conditions, but the turbulent boundary layer is then a 
relatively unimportant contributor to cabin noise. 
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1. INTRODUCTION 

High speed turbulent flow over an aircraft fuselage is responsible for a substantial 
component of the interior noise [1-4], and is probably the most important source of cabin 
noise for jet powered passenger aircraft in steady cruise. The principal mechanism of sound 
generation in subsonic flight is the interaction of boundary layer pressures with structural 
inhomogeneities, such as ring-stiffeners, stringers and windows, which produce sound whose 
amplitude and frequency content typically vary widely with flight speed, even though the 
boundary layer characteristics are believed to be only weakly dependent on flow velocity, at 
least within the subsonic regime. 

Several numerical schemes are being developed to predict the boundary layer generated 
interior noise [5]. Graham [6, 7] has proposed a semi-analytical model that is applicable to 
turbulent flow over an isolated fuselage panel, and which incorporates important structural 
effects of interior trim. None of these methods explicitly includes the influence of the 
exterior mean flow in calculating the noise. A simplified treatment that involves mean flow 
was proposed by Howe and Shah [8], who considered the interaction of a boundary layer 
with a nominally smooth elastic plate with periodically spaced rib-stiffeners. An elementary 
model of “coincidence scattering”, i.e., of sound generation by the edge scattering of flow 
excited structural waves whose phase velocity exactly matches the convection velocity of a 
flow inhomogeneity, was discussed in [9], but no account was taken of mean flow. However, 
although there is currently no fully satisfactory prediction scheme, the basic principles 
appear to be well understood. 

Most numerical prediction procedures take some account of lateral curvature of aircraft 
panels (e.g., see [5, 10]) and of interior panel trim, but there is no quantified understanding 
of the influence on sound generation of the step-like discontinuities in the outer skin that 
occur where adjacent fuselage panels overlap. This might be important when the step 
height is comparable to the boundary layer displacement thickness, so that the interaction 
between the turbulence and step is likely to be strong. Sound is generated by the scattering 
of turbulence pressures by a step, but in addition, the boundary layer wall pressures 
downstream of the step are also modified, and this will in turn affect the production of 
sound by direct panel forcing (e.g., see [11, 12]). 

In this chapter an analysis is performed to estimate the relative contribution from skin 
steps to boundary layer interior noise. The method is a development of that used by Howe 
and Shah [8] for a ribbed panel, which explicitly incorporates the influence of the finite 
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Mach number of the boundary layer. No account is taken in this calculation of the changes 
in boundary layer structure caused by the step, which is therefore assumed to generate 
sound by interaction with a pre-existing turbulent field that is swept over the step, as 
in the analogous problem discussed in [13] for a step in a rigid wall. Such changes are 
associated with separation induced at the step and will be important when the step height 
is comparable to the boundary layer displacement thickness. From an acoustic point of 
view separation tends to reduce the radiated sound relative to that predicted when the 
unsteady motion is modeled according to the proposed inviscid scattering theory [14, 15], 
because “noisy” potential flow singularities are smoothed-out. The principal conclusion 
of this investigation is that, for typical aircraft structures and cruise Mach numbers, the 
influence of skin steps on interior noise is small except at very high frequencies, which are 
probably too large to be relevant in practice. 

The analytical model is formulated in §2, where the sound generated by the interaction 
of the turbulent flow with a skin joint is expressed in terms of the unsteady, boundary layer 
wall pressure and a suitable acoustic Green’s function that includes the combined effects of 
structural compliance, the skin step and a rib stiffener at the panel overlap. This is used 
(in §§3,4) to investigate the dependence on Mach number of skin step noise relative to the 
noise generated by interaction of the boundary layer with the clamp line of adjacent panels. 
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2. THE AERODYNAMIC SOUND PROBLEM 


2.1 Formulation 

Consider high subsonic mean flow over the idealized skin step depicted schematically 
in Figure 1(a), which consists of two skin sections, modeled as semi-infinite, thin elastic 
plates which overlap and are clamped together along a nominally rigid, rectilinear rib or 
stringer. The exterior fluid has mean density pi and sound speed c x with mean flow in the 
aq-direction of the rectangular coordinate system (xi,x 2 ,x-z), as indicated in the figure. 
The overlapping skin sections are in practice riveted along the clamp line, which is taken to 
coincide with the 13-axis. 

When the influence of the skin step on sound production is ignored, the skin thickness h, 
say, is taken to be sufficiently small that both plates may be assumed to lie in the plane 
x 2 — 0, where x 2 > 0 in the exterior mean flow. This flow is assumed to be turbulent, and 
the problem of calculating the sound produced by the interaction of the turbulence with 
the plates in the neighborhood of the clamp line was considered in [8]. The basis of this 
method will be briefly recalled, and the extension necessary to take account of the skin step 
will then be described. 


The fluid in the “interior” domain x 2 < 0 has mean density and sound speed respectively 
equal to p 0 and c 0 and is in a mean state of rest. To determine the sound radiated into 
this region, the configuration of Figure 1(a) is first taken in the simplified form illustrated 
in Figure 1(b), where the skin step is ignored, and the adjacent skin sections are regarded 
as smoothly clamped together along the a; 3 -axis (at O). If U denotes the mean flow speed 
(in the aq-direction) in the exterior region outside the turbulent boundary layer, where the 
Mach number M = U(c\ < 1, the equations describing the production of aerodynamic 


sound can be taken in the form 


A & 2 


- V‘ 


B = 0, x 2 < 0, 


( 1 ) 


1 


r 2 

C 1 



B 


Q(x, t), x 2 > 0, 


( 2 ) 


where B = w + \v 2 is the total enthalpy, w being the fluid specific enthalpy, v the velocity, 
and t denotes time [16]. The source term Q(x, t) is non-zero only within the boundary 
layer, and includes all effects associated with the boundary layer mean shear flow and the 
turbulence sources. 
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When the motion is isentropic, the acoustic pressure fluctuation p can be obtained from 
the solution of equations (1), (2) by means of the relation dp/dt = pDB/Dt, where p is the 
local mean density and D/Dt is the material derivative. In those regions where the flow is 
irrotational, B = —d<t>/dt, where (j>(x,t) is the velocity potential of the unsteady motion. 

Deflections of the plates are assumed to be governed by the linearized thin plate equation 


D 


d 2 d 2 \ 2 d 2 

+ +m d? 


C+M = o, 


(3) 


where £(21, £3,2) is the flexural displacement in the a^-direction, D the plate bending 
stiffness, m the mass per unit area, and 


[p] = p{x u + 0 ,x 3 ,t) -p(x i, — 0, x 3 , t) 


is the pressure loading. The displacement C and the pressure are also related by the 
linearized kinematic relations 


d 2 ( _ —1 dp 
dt 2 p 0 dx 2 ' 


x 2 = -0; 



x 2 — +0. 


Along the clamp line the displacement satisfies the clamped edge conditions 


(4) 


c = 0, = 0, x 3 = 0. (5) 

OX i 

2.2 Green’s function in the absence of the skin step 

Green’s function G 0 (x, y,t — r) for the coupled equations (1), (2) is the solution with 
outgoing wave behavior, which also satisfies the appropriate conditions on the plates and on 
the clamp line, when Q on the right of (2) is replaced by the point source <5(x — y)S(t — r) 
concentrated at y. When G is known the acoustic pressure in the interior region (where 
B s; p/po) can be written 

p(x, t) = PoJ G 0 (x, y, t - t)Q( y, r)d 3 ydT, (6) 

where the spatial integration is over the boundary layer source region and the time 
integration is over — oo < r < oo. 
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When the motion is isentropic, the acoustic pressure fluctuation p can be obtained from 
the solution of equations (1), (2) by means of the relation dp/dt — pDB/Dt, where p is the 
local mean density and D/Dt is the material derivative. In those regions where the flow is 
irrotational, B = —d<j)/dt, where <£(x, t) is the velocity potential of the unsteady motion. 

Deflections of the plates are assumed to be governed by the linearized thin plate equation 


D 


( d 2 

d 2S 

\ 2 a 2 ' 

\foi 

+ dxl , 

+ m w_ 


C + [p] = 0, 


( 3 ) 


where ( \{x\,xz,t ) is the flexural displacement in the ^-direction, D the plate bending 
stiffness, m the mass per unit area, and 


[p] = p{x i , +0, rr 3 , t) - p{xi , -0, z 3 , t) 


is the pressure loading. The displacement £ and the pressure are also related by the 
linearized kinematic relations 


d 2 C, — 1 dp 

dt 2 p 0 dx 2 



x 2 — +0. 


Along the clamp line the displacement satisfies the clamped edge conditions 


( 4 ) 


C = 0, — = 0, :r 3 = 0. (5) 

OX\ 

2.2 Green’s function in the absence of the skin step 

Green’s function G 0 (x,y,t — r) for the coupled equations (1), (2) is the solution with 
outgoing wave behavior, which also satisfies the appropriate conditions on the plates and on 
the clamp line, when Q on the right of (2) is replaced by the point source <5(x — y )S(t — r) 
concentrated at y. When G is known the acoustic pressure in the interior region (where 
B « p/Po) can be written 

P(x, t) = p 0 J G 0 (x, y, t - t)Q( y, r)d 3 ydT, (6) 

where the spatial integration is over the boundary layer source region and the time 
integration is over — oo < r < oo. 
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It will suffice for the present discussion to confine attention to the generation of sound 
by boundary layer disturbances Q = Q(x u x 2 , t) that are independent of the spanwise 
coordinate x 3 . This is a special case of the fully three-dimensional analysis given in [8], and 
is adequate for estimating the relative contribution to the sound from the skin step. For an 
observer at x in the interior domain, and at large distances from the clamp line, we then 
have (see [8] for details) 


C 0 (x,y, t — t ) 


Vipx/Po f°° T( a;, (9) ( - itio{niyi-N 2 y2 ) 
27t(87t|x|) 1 / 2 Joo V 


+ 


No 


2, ( a> + C/fe)o n A:"e»( fc » I+r+(fc)lB) dk ^-i^t-r 


£/- 


cxico/d - M sin 9) £ \J-oo T + (k)C + (k,u) 


-W c ° ] du, 
|x| — > oo, 


( 7 ) 


where the integration contour in the fc-plane is indented to pass respectively above and 
below real poles and branch points of the integrand in k > 0, T(u>,9) is the elastic plate, 
plane wave transmission coefficient 


T(lo,9) 


( j* — M sin 



( 8 ) 


and the remaining new terms appearing in these formulae are defined as follows: 


m = sin6>, N 2 = y/(%~ Msin^) 2 - sin 2 0, e = p 0 c 0 /u c m, u c = c 2 0 (m/D ) 1/2 , 

9 being the angle shown in Figure 1(b) defining the observer direction, 

7 (jfe) = (k 2 0 - k 2 ) 1/2 , T ± {k) = [(c j/c, ± Mk ) 2 - fc 2 ] 1/2 , /c 0 = u/c 0 , 

£±(k,uj) = Dk* - mu 2 - {[poL} 2 /*/^) + pi(w ±Uk) 2 /r±(k)} , (9) 

where branch cuts in the upper and lower halves of the complex A;-plane are taken such 
that 7 (lt) and r ± (A;) have positive imaginary parts on the real axis. Here, u c > 0 is the 
coincidence frequency above which the in vacuo phase speed of flexural waves on the plate 
exceeds the interior sound speed c 0 . The coefficient e is a measure of the influence of fluid 
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loading on the motions of the plates. For an aluminum plate of thickness h in air, e « 0.0022 
and u c h/c 0 « 0.22. The four coefficients a n are complex functions of the frequency w, whose 
values depend on conditions at the junction of two plates, and satisfy a n (—u>) = a* (a;), 
where the asterisk denotes complex conjugate. For the clamped condition considered here 
[ 8 ] 


do 


mu 

~kI 


Aq, Oi 


muj 2 

~Kf 


A\, a 2 = 0, a 3 = 0, 


( 10 ) 


where K 0 = (mcj 2 /D) 1/4 > 0 is the in vacuo wavenumber of bending waves of frequency u>, 
and the dimensionless coefficients Ao, A\ are given by 


^2 + fJ-n 1^1 


_ — (^i + /zni^p) 

1 _ - 'I'o^ 



(11) 




w: 


X j dX 


A 4 - 1 - 




+ 


(Pocg/Pic?)(co/ci+ MA /M) 2 A 

y/UtCo/a+MX)*-* 1 ) . 


( 12 ) 


The path of integration along the real axis is indented to pass respectively above and below 
poles and branch points of the integrand in A £ 0. 


2.3 Skin step modified Green’s function 

In the presence of the skin step the representation (6) becomes 


p(x, t) = p 0 J |G 0 (x,y,t- t) + G s (x,y,t-r)J(2(y > T)d 3 ydr, (13) 


where G s is the correction to Green’s function (7) produced by the presence of the step. To 
calculate G s we first recall that the reciprocal theorem [8, 17, 18] implies that the term in 
the large brackets of the integrand of (7), regarded as a function of y, is proportional to the 
velocity potential in y 2 > 0 generated by a point source at the observer location x when 
the mean flow in the exterior region is reversed in direction. This potential must satisfy 
boundary conditions of the type (3) - (5) in the reverse flow problem. Thus, if we assume 
the step to be in the immediate neighborhood of the clamp line 2/3 = 0 (as in Figure 2), the 
presence of the step implies that the normal surface derivative dG/dy n = d(G 0 + G s )/dy n 
should vanish on the step profile in the reversed flow problem. In the reciprocal problem 
the length scale of the unsteady fluid motion in the neighborhood of the step when the 
step is ignored is of order l/K 0 < 1/k 0 . If K 0 h < 1, i.e., if the structural wavelength is 
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large compared to the skin thickness, which is usually the case (and is actually necessary to 
justify the use of the thin plate equation (3)), the velocity potential of the unsteady motion 
in the neighborhood of the step relative to that of the plates is equal to a sum of terms 
proportional respectively to ni[y L + <£*(y)] and k[yi + <£*(y)], where <f>*( y) is a solution 
of the homogeneous, time harmonic form of equation (2) when the mean flow is reversed , 
which satisfies 


^ = 0 on the step profile. (14) 

9]Jn 

This is because, near the step, the exponentials in (7) can be expanded to first order in y; 
the term in yi must then be augmented by <f>* in order to satisfy (14). If the coordinate 
origin is taken at the foot of the step, then to first order in the step height s (which is equal 
to plate thickness h for simple overlapping plates), the boundary condition (14) becomes 


df_ 
dy 2 


= ±S(5(yi) on y 2 = +0, 


and the appropriate <f>* is accordingly found to be given by 

±s r°° e*(* !,1+r+ (*) y2 ) 


0*(y) = — f 

* ' 2 t n J-c 


r +(*) 


dk, 


(15) 


where the ± sign is taken according as the step is forward or backward facing relative to 
the mean flow. We then have 


G s (x,y,t - r) 


(Pi/po)<ft*(y) r 00 

2nsJ&Tci\x\ 


Tju^e) 


{ 


K 0 n i 


No 


Ci(c 0 /ci - MsmO) ^ 0 I-oo r + (k)£+(k,u>) ] 


(q> + f/fc)q n /c n+1 dk 1 -ia;(t-r-|x|/co) 


dui. 


(16) 


The integral within the brace brackets converges only for n < 1; in particular it converges 
for the clamped condition, for which a 2 = 03 = 0. This is because when a 2 7^ 0 there is a 
discontinuity in the slope of the two plates at the join, and there can then be no smooth flow 
over the plates in the absence of the step (see [8]), as assumed in the reciprocal problem. 
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3. THE INTERIOR ACOUSTIC PRESSURE 

3.1 The blocked wall pressure 

It is convenient to express the interior radiation in terms of the boundary layer wall 
pressure on the exterior skin that would be generated when the wall is regarded as rigid. 
This is the blocked pressure, and can be calculated by solving (2) subject to the boundary 
condition dB/dx 2 = 0 on x 2 = +0, and using the relation dp/dt « pi(d/dt + Ud/dxi)B. If 
is the component of the blocked pressure of frequency u> and wavenumber 
ki in the streamwise direction, then it is readily shown that 


where 


Pb(k\,u) 


2nipi(u> — Uk\) 
uT-iki) 


Q(^ii — r_ (fci ), oj) , 


(17) 


Q(k u k 2 ,uj) = ^^3 JJ dx\dt Q(x u x 2 ,t)e l{klXl+k2X2 ut) dx 2 . (18) 

The dominant components of the blocked wall pressure occur in the vicinity of ki = u>/U c 
where U c is a convection velocity approximately equal to 0.7 U [19 - 22]. 

3.2 Sound radiated from the clamp line 

The acoustic pressure p 0 (x, t), say, produced by the interaction of the boundary layer 
pressure field with the clamp line is given by (6). The contribution from the first, 
exponential term in the round brackets in the representation (7) of G 0 can be ignored, since 
it does not correspond to the generation of sound at the clamp line, but to the transmission 
through the wall of boundary layer pressure fluctuations that are already sound, exactly as 
if the clamp were absent. Relations (17) and (18) can be used to express the contribution 
from the second term in the round brackets of (7) in the form 


PoCM) 


71 


N 2 


2i|x| CL (^-M sine) 

Co \ Cl ) 


t r r(a; ’ e y^ a ^ - ^ n - p b (k u cje-MHxi/*) dkldU ' 

n^Q J-oo £-{k i,(Jj 


(19) 


3.3 Sound radiated from the skin step 

Similarly, the contribution p s (x,t) to the acoustic pressure from the step, given by the 
term in G s in (13) becomes 
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p.(x,t) * 5£./4l / /" 

27r Y 2i|x| J 7-oo (u — Uki) \ 
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y> f°° + Uk)ank r, ~ il dh \ i -iu(t-|x|/e„) j.j 

„V-» r + (*)£ + (*, u ) P 1 * 1 ’"'' 


KoCi — A/ sin$) n=0 J-oo r + (/c)£ + (fc 


( 20 ) 


the upper/lower sign being taken respectively for forward and backward facing steps. 

3.4 Acoustic efficiencies 

Consider the sound generated by the interaction of the component pi = pb(ki,u)dkiduj 
of the blocked wall pressure with the clamp line and skin step. The relative efficiencies of 
sound generation by these two mechanisms is determined by the relative magnitudes of the 
ratios p 0 (x.,t)/pi and p s (x, t)/pj. Let us suppose that pi is characteristic of the dominant 
wall pressure fluctuations, for which u/k i = U c . It then follows from (19) and (20), using 
also results given in §2, that 
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4. NUMERICAL RESULTS 

The curves in figures 3 - 5 are plots of the calculated efficiency 20*log (|p(x, t)/pj\/,JhJ\xf) 
(dB) of interior sound production by the clamp line and by the skin step at three different 
mean flow Mach numbers M, and for the nominal radiation direction 9 = 45° (see Figure 
1(b)). The calculations have been performed for aluminum plates in air; for simplicity it 
is assumed that p 0 = Pi, c 0 = Ci, for which the remaining non-dimensional parameters 
describing the fluid-structure interaction have been assigned the values e = 0.0022, 
u) c h]c 0 — 0.22, M c = 0.7 M. 

The peaks labeled “convective resonance” occur where the convection velocity u/ki — U c 
of the blocked wall pressure component Pb{ki,cj) is the same as the phase speed of bending 
waves in the plates of the same frequency and wavenumber, and occurs approximately at 
ui = M^u> c . At frequencies exceeding the convective resonance frequency the efficiency of 
sound generation by scattering at the clamp line decreases very rapidly. At frequencies 
below resonance the overall level of the skin step generated sound is always negligible. At 
higher frequencies, however, it becomes dominant because of the precipitous decrease in the 
efficiency of sound generation by clamp line scattering. The peaks in the skin step sound 
evident in Figures 4 and 5 occur at sufficiently large frequencies that u)fu) c = 1/sin 2 9 > 1. 
This condition is satisfied by those structural waves generated at the step whose phase 
velocity is just equal to the trace velocity of sound waves launched in the observer direction 
6, and corresponds to the vanishing of the first term in the brace brackets of the denominator 
of the definition (8) of the plane wave transmission coefficient T(u,6). 

The figures indicate that the skin step sound will be dominant at sufficiently high 
frequency. To illustrate the orders of magnitude involved, Table 1 lists approximate values 
of the critical frequency uj s /2-k = f a Hz, say, above which skin step noise is dominant for the 
three Mach numbers considered in Figures 3-5 when 6 = 45°, and when the skin thickness 
h = 0.1 cm. 
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Mach Number 

uj s /u) c 

fs Hz 

0.01 

0.0008 

10 

0.4 

0.17 

2000 

0.8 

1 

12000 


Table 1 . Critical frequency when 6 = 45°. 


Thus, at M ~ 0.8, which is appropriate for cruise conditions at which boundary layer 
generated interior noise is believed to be significant, it appears that the contribution from 
skin steps will be unimportant at frequencies less than, say, 10 kHz. According to Table 1, 
the critical frequency decreases rapidly with decreasing flight speed, being as small as 2 kHz 
when M = 0.4. However, such flight speeds are characteristic of takeoff and landing, when 
boundary layer generated noise within the cabin is negligible compared to that generated by 
other sources. The critical frequency varies relatively little with radiation direction 6. This 
is evident from Figure 6, which depicts the variation of 20 x log \p 0 /p s \ (dB) with frequency 
for 9 = -45°, 0°, 45°. 
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5. CONCLUSION 

Turbulent boundary layer generated noise is one of the dominant components of passenger 
aircraft interior noise at high subsonic cruise conditions. This noise is produced by 
the interaction of turbulent pressures with fuselage inhomogeneities, in particular with 
skin panel junctions. In this chapter an analysis has been made to estimate the likely 
contribution from interaction of the turbulent flow with skin steps, where adjacent panels 
overlap. Sound and structural vibrations are excited by the interaction. In addition, the 
step causes a localized, temporary change in the boundary layer turbulence characteristics. 
This has been ignored, and the turbulence has been regarded as frozen during its interaction 
with the step. Skin step noise is found to be significant only at high frequencies, beyond the 
structural convective resonance frequency (at which the boundary layer eddy convection 
velocity coincides with the phase speed of flexural waves in the skin); above this frequency 
the efficiency of sound production by the scattering of turbulence pressures at a clamped 
panel edge decreases precipitously. However, at cruise Mach numbers, where boundary layer 
generated noise is important, the skin step noise is dominant only at frequencies typically 
in excess of about 10 kHz, which is beyond the range of interest in practice. At lower flight 
speeds, appropriate to take-off and landing approach conditions, the skin step contribution 
is significant at lower frequencies, but in these circumstances cabin interior noise is no 
longer controlled by the fuselage boundary layer. 


75 


Report No. AM-97-027 


Boston University, College of Engineering 


REFERENCES 

1. J. F. Wilby, C. D. McDaniel and E. G. Wilby 1985 NASA CR-178004. In-flight acoustic 
measurements on a light twin-engined turboprop airplane. 

2. J. S. Mixson and J. F. Wilby 1991 in Aeroacoustics of Flight Vehicles: Theory and 
Practice, Vol. 2. NASA Ref. Pub. No. 1258, Chapter 16: Interior noise. 

3. W. R. Graham 1993 Boundary layer noise and vibration. Doctoral Thesis, Cambridge 
University Engineering Department. 

4. Wilby, J.F. 1996 Journal of Sound and Vibration 190, 545 - 564. Aircraft interior 
noise. 

5. R. J. Silcox 1995 (editor) Proceedings of the Interior Noise Workshop, held at NASA 
Langley Research Center, 25 - 27 April. 

6. W. R. Graham 1996 Journal of Sound and Vibration 192, 101 - 120. Boundary layer 
induced noise in aircraft, Part I: the flat plate model. 

7. W. R. Graham 1996 Journal of Sound and Vibration 192, 121 - 138. Boundary layer 
induced noise in aircraft, Part II: the trimmed flat plate model. 

8. M. S. Howe and P. L. Shah 1996 Journal of the Acoustical Society of America 99, 3401 
- 3411. Influence of mean flow on boundary layer generated interior noise. 

9. P. L. Shah and M. S. Howe 1996 Journal of Sound and Vibration 197, 103 - 115. Sound 
generated by a vortex interacting with a rib-stiffened elastic plate. 

10. W. R. Graham 1995 Journal of the Acoustical Society of America 98, 1581 - 1595. The 
influence of curvature on the sound radiated by vibrating panels. 

11. T. M. Farabee and M. J. Casarella 1988 ASME paper entitled: Wall pressure fluctua- 
tions beneath a disturbed turbulent boundary layer, pp 121 - 135 of Acoustic Phenom- 
ena and Interaction in Shear Flows Over Compliant and Vibrating Surfaces-Volume 6 
(Editors: W. L. Keith, E. M. Uram and A. J. Kalinowski.) 

12. J. Sulc, J. Hofr and L. Benda 1982 Journal of Sound and Vibration 84 , 105 - 120. 
Exterior noise on the fuselage of light propeller driven aircraft in flight. 

13. M. S. Howe 1989 Journal of Fluids and Structures 83 , 83 - 96. Sound produced by 
turbulent boundary layer flow over a finite region of wall roughness, and over a forward 
facing step. 


76 


Report No. AM-97-027 


Boston University, College of Engineering 


14. M. S. Howe 1976 Journal of Fluid Mechanics 76, 711 - 740. The influence of vortex 
shedding on the generation of sound by convected turbulence. 

15. M. S. Howe 1988 Proceedings of the Royal Society A420, 157 -182. Contributions 
to the theory of sound production by vortex-airfoil interaction, with application to 
vortices with finite axial velocity defect. 

16. M. S. Howe 1975 Journal of Fluid Mechanics 71, 625 - 673. Contributions to the 
theory of aerodynamic sound, with application to excess jet noise and the theory of 
the flute. 

17. M. S. Howe 1975 Journal of Fluid Mechanics 67, 579 - 610. The generation of sound 
by aerodynamic sources in an inhomogeneous steady flow. 

18. W. Mohring 1979 Modelling low Mach number noise, pp 85 - 96, Proceedings of the 
Symposium on Mechanics of Sound Generation in Flows , Gottingen, August 28 - 31 
(editor E.-A. Muller: Springer, Berlin). 

19. B. M. Efimtsov 1982 Soviet Physics Acoustics 28 , 289 - 292. Characteristics of the 
field of turbulent wall pressure fluctuations at large Reynolds numbers. 

20. D. M. Chase 1987 Journal of Sound and Vibration 112 , 125 - 147. The character of 
the turbulent wall pressure spectrum at subconvective wavenumbers and a suggested 
comprehensive model. 

21. A. E. Landmann 1995. Spatial-Temporal boundary layer models for aircraft interior 
noise. Paper presented at the NASA Langley Workshop on Interior Noise, Hampton 
VA, 25 - 27 April. 

22. W. R. Graham 1995 AIAA Paper 95-097. A comparison of models for the wavenumber- 
frequency spectrum of turbulent boundary layer pressures. 


77 



Report No. AM-97-027 


Boston University, College of Engineering 


(a) 


*2 

A 

x 3 


u 


exposed skin step 

/ 


rib/stringer 


(b) 



Figure 1. (a) Schematic skin-step configuration. 

(b) Simplified model in absence of skin step. 
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Figure 3. Efficiencies of sound generation by scattering from the clamp line and from 
the skin step; aluminum skin, M = 0.01, 0 = 45°. 
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Figure 4. Efficiencies of sound generation by scattering from the clamp line and from 
the skin step; aluminum skin, M = 0.4, 6 = 45°. 
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Figure 5. Efficiencies of sound generation by scattering from the clamp line and from 
the skin step; aluminum skin, M = 0.8, 6 = 45°. 
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Figure 6. Ratio |p 0 /p s | 2 (in dB) of the sound generated respectively by scattering 
at the clamp line and skin step for different radiation angles 9\ aluminum 
skin, M — 0.8. 
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CHAPTER 4 

INFLUENCE OF SEPARATION ON SOUND 
GENERATED BY VORTEX-STEP INTERACTION 
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SUMMARY 

An analysis is made of the sound produced when a line vortex interacts at low Mach 
number with forward or backward facing steps. The radiation is dominated by an 
aeroacoustic dipole whose strength is equal to the unsteady drag on the step. The drag is 
determined by the vorticity distribution, and a correct estimate of the sound must therefore 
include contributions from vorticity in the separated flow induced by the vortex. Separation 
is modelled by assuming that shed vorticity rolls up into a concentrated core, fed by a 
connecting sheet from the edge of the step of negligible circulation. The motion everywhere 
is irrotational except at the impinging vortex and the separation core, and the trajectory 
of the core is governed by an emended Brown and Michael equation. For large steps it is 
found that estimates of the generated sound that neglect separation are typically an order of 
magnitude too large. The sound levels predicted for small steps with and without separation 
are of comparable magnitudes, although the respective phases are different. Turbulent 
flow over a step frequently involves separation and large surface pressure fluctuations 
at reattachment zones. The results of this chapter suggest that numerical schemes for 
determining the noise generated by turbulent flow over a step must take proper account of 
“forcing” of the separation region by the impinging turbulence and of vorticity production 
via the no-slip condition. 
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1. INTRODUCTION 

Boundary layer turbulence and other flow inhomogeneities can induce unsteady separation 
from structural irregularities, such as rivet heads, lap-joints, cutouts, etc. The unsteady 
drag exerted on the irregularity is equivalent to a localized aeroacoustic dipole, and is 
responsible for the production of both sound and structural vibration. Lap-joints (or “skin 
steps”) occur at junctions of neighboring panels on an aircraft fuselage, and are a possible 
source of boundary layer generated interior cabin noise, either because they enhance 
downstream wall pressure fluctuations that drive vibrations of the cabin walls, or because 
the drag dipole is a direct source of sound and vibration (Efimtsov 1996; Howe 1997). 

Irregularities in the form of forward or backward facing steps have been studied extensively 
(Bradshaw and Wong 1972; Chandrsuda and Bradshaw 1981; Eaton and Johnston 1981; 
Farabee and Casarella 1984, 1986, 1988; Moss and Baker 1980) because their free shear 
layers possess well defined separation and reattachment zones. More generally, reattachment 
occurs downstream of a stagnation point, and is a region of high fluctuating wall pressure, 
where free shear layer eddies bifurcate into upstream and downstream travelling components 
(Bradshaw and Wong 1972). According to Farabee and Casarella (1984, 1986), low Mach 
number turbulent wall pressures at reattachment for backward and forward facing steps 
(with step heights comparable to the boundary layer thickness) are respectively about 5 
and 10 times larger than the smooth wall boundary layer pressures. The wall pressure 
must ultimately revert to that of a smooth wall boundary layer, but differences have been 
observed at distances as large as 72 step heights downstream of the step. 

The large changes in the boundary layer structure at a step must be carefully modeled 
to obtain accurate estimates of the sound produced by the boundary layer-step interaction. 
For example, in low to moderate Mach number flows, the unsteady drag (which is the 
dipole source strength) is determined by the vorticity distribution near the step (Howe 
1989a, 1995), which includes vorticity in the impinging flow together with any produced by 
the interaction and “trapped” in the recirculating separation bubble. Hitherto, theoretical 
predictions of the generated sound (e.g., Conlisk and Veley 1985; Dhanak and Gundlapalli 
1992; Howe 1989b, 1997) have tended to ignore separation and vorticity production. Conlisk 
& Veley (1985) and Dhanak & Gunlapalli (1992) represented the impinging inhomogeneous 
flow by an assembly of line (or point) vortices, whose motions over a large step (comparable 
in size to the boundary layer thickness) were determined from inviscid equations of motion. 
Howe (1989b, 1997) considered “small” steps, and calculated the sound as a by-product of 
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the scattering of the boundary layer pressure field, which was assumed to be frozen during 
convection over the step. 

In this chapter we assess the accuracy of these inviscid approximations from an analysis of 
the interaction of a single line vortex with a forward or backward facing step. Separation is 
modeled by assuming that vorticity shed during the interaction rolls up into a concentrated 
core, which grows in strength by the passage of continuously shed vorticity along a 
connecting sheet from the edge of the step. The motion is irrotational everywhere except 
at the impinging vortex and the core of the shed vorticity, and the trajectory of the core 
is determined by the emended Brown and Michael equation (Brown and Michael 1954, 
1955; Howe 1996). For large steps, it is concluded that estimates of the sound that neglect 
separation are typically an order of magnitude too large. This is because the calculated 
drag on the step is greatly overestimated when the potential theory singularity at the edge 
of the step is not suppressed by vortex shedding. For small steps, predicted sound levels 
with and without separation are of comparable magnitudes, but exhibit large differences in 
phase. 

The aerodynamic sound problem is formulated in §2 for low Mach, self-induced vortex 
motion over forward or backward facing steps on a plane wall. Numerical predictions of the 
sound are presented in §3. 
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2. THE AERODYNAMIC SOUND PROBLEM 

2.1 Formulation 

A line vortex of circulation T 0 > 0 is in translational motion adjacent to a rigid wall. The 
vortex is parallel to the z-axis of the rectangular coordinate system ( x,y,z ), and the wall 
is parallel to the plane y = 0, except for a “vertical” step of height h, whose foot coincides 
with the z-axis. To fix ideas, consider the forward facing step illustrated schematically in 
Figure 1, where the coordinate origin is at O, and the fluid and vortex lie in the region 
y > 0 “above” the wall. 

In the undisturbed state the fluid is at rest, and is inviscid except that vortex shedding 
is permitted from the right-angled edge S of the step to remove the singular velocity and 
pressure that would otherwise occur. Vortex shedding is uniform along the step, and is 
modeled by a line vortex of variable strength T(t) ( t denoting time) whose axis intersects the 
plane z = 0 at ( x , y) = x r (f), and translates at velocity dx r /dt = V. The vortex strength 
increases as shed vorticity is continuously “fed” along a connecting sheet from S. The 
circulation of the connecting sheet is assumed to be negligible compared to T. This model 
of vortex shedding, involving a time-dependent core T containing all of the shed vorticity, 
includes a pressure jump across the connecting sheet, and requires the translational velocity 
V to be different from the fluid velocity v at the core (Brown & Michael 1954, 1955). T(f) 
is assumed to vary monotonically until the vortex is finally “released” from S when dT/dt 
changes sign. The “free” vortex (with no connecting sheet) then proceeds to convect with 
the fluid at velocity v, and a new vortex is released from S. In the present discussion dT/dt 
turns out to be one-signed during the dominant interaction of the incident vortex r o with 
the step. 

Let (x 0 , y Q ) be the position of the incident vortex in the plane z = 0 at time t. The 
vorticity w can then be written 

u) — D 0 4* fl, 

D 0 = r o k5(a:-T 0 (t))(5(y-yo(f)), 

D = T(t)k8{x-x r (t))5(y-y r (t)), (2.1) 

where k = (0, 0, 1) is a unit vector in the z-direction (out of the plane of the paper in Figure 
1). For an ideal fluid the momentum equation is 
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l^(/f + H = -“ Av + 7 (2 ' 2) 

where p and p are respectively the pressure and density, and for two-dimensional motion 
v = (u, v, 0). F is the distributed force 


F = pfl A (v — V) — p— n5(xj_)H(s r - s), s > 0, (2.3) 

where x± is distance from the connecting sheet measured in the direction of the local normal 
n of Figure 1, s r is the length of the sheet (between S and T), and s is distance measured 
along the sheet from S. The first term on the right is concentrated at the core of the shed 
vortex, and is equal and opposite to the Joukowski lift on T which arises when V^v; the 
second term represents the pressure force across the connecting sheet (Howe 1996). 

2.2 Equations of motion of the vortices 

The incident vortex is “free” and translates at the fluid velocity v(x, y , t) evaluated at its 
core 


dx o 

dt 


= v(x 0 ,y 0 ,t). 


(2.4) 


This velocity is calculated according to potential flow theory, and consists of a contribution 
from the images of r o in the rigid wall and from the shed vortex T and its images. 


Brown and Michael (1954, 1955) obtain the equation of motion of the shed vortex by 
requiring the net force exerted on the fluid by the distribution F to vanish. In general, 
it is not sufficient to require / F dxdy - 0, since the surface reaction to the remaining 
unbalanced couple (/ x A F dxdy ) is usually equivalent to a surface force on the fluid. 
However, for an infinitely large rigid wall, the aggregate normal reaction is equal and 
opposite to the normal (i.e., the y-) component of /F dxdy , and it is therefore sufficient to 
require that this should vanish. Performing the integration for the y-component of (2.3), we 
find 


dx r 
dt 


x r dT , .. 

+ Y~dt = u ( Xr ’3 /r ’ f )‘ 


(2.5) 


The net horizontal component of the force on the fluid (in the ^-direction) induced by F 
and its interaction with the step is eliminated by requiring the shed vortex path to satisfy 
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the equation (Howe 1996) 

j F VX dxdy = 0, (2.6) 

where X = X(x,y) is a harmonic function that depends only on the shape of the step, and 
is equal to the velocity potential of ideal, incompressible flow past the step that has unit 
speed in the x-direction at large distances from the step. Equation (2.6) is expressed in 
differential form by making use of the Cauchy-Riemann relations 

dX/dx = dV/dy, dX/dy = -d^/dx, (2.7) 

where is the stream function conjugate to X. When ^ = 0 on the wall, equation (2.6) is 
equivalent to 


dxr „ 'I' dr . . . , /n 

___ . + -— = v • V^, (x,y) = (x r ,y r )- (2.8) 

This equation is reduced to the following simplified form by subtracting equation (2.5) 
multiplied by d^/dx 


dyr (fr - x r fr x ) 1 dr 
dt T dt 


where (^ x , ^ y ) = (d^ /dx,d^ /dy). 


u(a:r,2/r,0> 


(2.9) 


2.3 The aerodynamic sound 

At low Mach numbers the production of aerodynamic sound by the vorticity and force 
distribution on the right of (2.2) is governed by the inhomogeneous wave equation (Howe 
1975, 1996) 


{d 2 j(? 0 dt 2 - V 2 )£ = div(u> A v) — div(F/p 0 ) (2.10) 

where B = p/p 0 + ^v 2 (the total enthalpy in isentropic flow) and p 0 , c 0 are respectively the 
mean density and sound speed. In irrotational regions the motion can be described by a 
velocity potential ^(x, t), and B = —d<f>/dt. Thus, at large distances from the step, where 
the perturbed motion is small, the pressure p « p 0 B, and the solution of (2.10) can be 
expressed in the form 

c BC 1 

p(x,t) « -J (p 0 w A V - F)(y,r) • — (x,y,t- r) d 3 ydr, y = (x',y',z'), (2.11) 
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where the integration is over the fluid and all times r. In this formula G(x, y,t-r) is the 
Green’s function with outgoing wave behavior that has vanishing normal derivative on the 
wall, i.e., the solution of (2.10) when the right hand side is replaced by <S(x - y)S(t - r). 

The characteristic wavelength of the generated sound is much larger than the step 
height h when the Mach number is small, and G is then well represented by the compact 
approximation (Howe 1989b) 


6 {t - r - |x - Y|/c 0 ) S ( t - T - |x — Y|/c 0 ) 
G(x, y,t — t) « 4tt|x- Y| + 4tt|x - Y| 


oo, 


Y = (. X{x\y'),y\z '), Y = (X(x',y'), -y',z'). 


( 2 . 12 ) 


To use this formula in the general solution (2.11), it will be assumed that the effective 
source region lies within the finite span -\i < z < \i, where l is small compared to 
the dominant acoustic wavelengths. This approximation is not essential to the following 
discussion, but corresponds more closely to conditions in practical problems, where the 
radiation spreads three-dimensionally rather than cylindrically. When end effects (from 
z = ±|£) of the source distributions are ignored, the only non-trivial contributions to the 
radiation are supplied by the following approximation to G, obtained by taking the first 
term in the expansion of (2.12) in powers of the retarded source position, 


G(x, y, t - t) « X( ^ x C ° S( V (< - r - M/C.) (2.13) 

where the prime denotes differentiation with respect to time, and 6 = cos _1 (x/|x|) is the 
angle between the observer direction in the acoustic far field and the x— axis. 

When this expression is substituted into the general solution (2.11) it follows from (2.6) 
and the Cauchy-Riemann equations (2.7), that the contribution from F vanishes identically, 
and that the acoustic pressure can be written 


p(x, t) 


p 0 £ cos 6 d 
27tc 0 |x| dt 


(i>- •v*) fcW + ( r <‘ )v -v*) I 


x r(0j t— | x |/co ’ 


:| — > oo, (2.14) 


where the notation implies that the terms in the square braces are evaluated at the retarded 
time t — |x|/c 0 . 
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3. NUMERICAL RESULTS 

3.1 Vortex trajectories 

The vortex paths are determined by equation (2.4) for the incident vortex r o and 
equations (2.5), (2.9) for the shed vortex T. The irrotational velocity v is computed in the 
usual way, by introducing the dimensionless complex variable 


and transforming the fluid region above the step onto the upper half of the complex C-plane. 
This is accomplished by the mapping 


Z = \ {\/C 2 - 1 =F In (C + v/C 2 - l) } + £(! ± 1), Im C > 0, (3.1) 

where here and henceforth, upper/lower signs are to be taken respectively for the 
forward/backward facing step. The foot of the step O maps into the point C = ^Fl, and the 
top S maps into C = ±1. 

Let Co, Cr denote the respective images in the C-plane of the incident and shed vortices. 
At low Mach numbers, when the step is acoustically compact, the motion near the step may 
be regarded as incompressible, with complex potential 

w = ^ ( ln(< - Co) - ln(c - o) + ^ ( ln(C - Cr) - MC - <?)) . (3.2) 

where the asterisk denotes complex conjugate. The instantaneous value of T(t) is determined 
from the Kutta condition that dw/dz should be finite at S, which is satisfied provided 
dw/dC, — > 0 as C — t ±1 respectively for forward and backward facing steps. This yields 



which implies that T and F 0 always have opposite signs. 

The velocity v(rr 0 , y 0 , t) of the incident vortex, on the right of equation (2.4), is calculated 
from the velocity potential (3.2) by first excluding the free field self-potential In (Z — Z a ), 

where Z 0 — x Q /h + iy 0 /h. In the C-plane equation (2.4) is equivalent to 
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where 


1 - - w«l 

\( 1 

o * 
1 


T = - 


+ 


1 


1 


1 


Co* -Co 2(C 0 2 -1)V r c \Q - Cf Co -Cr, 


2 h? 


K ± I 

Cti' 


(3.4) 


(3.5) 


Similarly, the velocity components u(xp,yp,t ) and v(xp,yp,t) in the equations of motion 
(2.5) and (2.9) of T are evaluated from (3.2) by excluding ^ ln(Z — Zp), Zp = xp/h + iyp/h. 
The stream function = Im(/iCA), so that, in the £-plane, Cr satisfies 


dCr 

dT 


+ 


R^))( ! Im(fr) + ” R e(2 r )^(C r ))i^ 


= -*l^(Cr)| 


r. (cf - Cr + rn - 1)V (c? - « 



(3.6) 


where Re(Zp) on the left hand side is given in terms of Cr by (3.1). 

Equations (3.3) - (3.6) are solved numerically by assigning a large and negative initial 
value to Re Co and adjusting the corresponding value of ImCo to make the stand-off distance 
d, say, of the vortex T 0 from the wall equal to some prescribed initial value. The starting 
value of Cr is taken at some point close to the image C = ±1 of the step top S; the 
precise location depends on the integration step length, but does not critically affect the 
solution, since T is very small when the shed vortex is close to S, and any starting errors 
rapidly become negligible as the solution builds up with the approach of r o to the step. A 
fourth-order Runge-Kutta procedure was used; the equations form a robust system and no 
difficulty was experienced in achieving convergence. 


3.2 Acoustic pressure 

The acoustic pressure (2.14) may now be expressed in the form 


p(x, t ) « &Tt 3 p 0 U 2 M 


d \ 3 i cos 6 d 
h) | x | dT 



+ 


r CO 
r 0 



t-|x|/Co 


M = —, (3.7) 
C 0 


where d^p/dT is proportional to the fluid velocity at the core of the shed vortex V (which 
is the same as the vortex translation velocity V only when dr jdt = 0) and is equal to the 
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right hand side of (3.6) evaluated at the retarded position of T. The velocity U = r o /4nd is 
the initial translational velocity (which is parallel to the wall) of the incident vortex before 
its motion is affected by the step. The acoustic pressure is proportional to p 0 U 2 M , which 
is typical of an aeroacoustic source of dipole type; the dipole strength is just the unsteady 
drag exerted on the step during the interaction. 

3.3 The forward facing step 

Consider first the hypothetical case of vortex motion and acoustic radiation in the absence 
of vortex shedding (Conlisk & Veley 1985; Dhanak & Gundlapalli 1992). Only the incident 
vortex r 0 is present; its trajectory is governed by equation (3.4) with T = 0. The trajectory 
is depicted in Figure 2a for d/h — i.e., when the initial distance of the vortex from the 

wall is half the step height. The points labeled along the path indicate the vortex position 
at various non-dimensional times Ut/h, where time is reckoned from the instant at which 
the vortex passes the step (at x = 0). 

When Ut/h — » — oo (far from the step) the kinetic energy of the (incompressible) motion 
per unit length of the vortex is equal to (p 0 Tl/4n) \n(d/r 0 ), where r c is a length of the order 
of the radius of the vortex core (Lamb 1932). In the absence of shedding this energy is 
conserved during interaction with the step, so that the distance of the vortex from the wall 
ultimately returns to d after the interaction when Ut/h becomes large. Because of this, the 
ratio Ut/h is approximately equal to the horizontal distance of the vortex from the step 
measured in step heights. 

The non-dimensional acoustic pressure 


p(x, t) 

p 0 U 2 M(d/h) 3 (£/\xj) cos0 ’ 

is plotted in Figure 2b as a function of the nondimensional retarded time U[t]/h, 

[f] = t — | x | / c 0 . Most of the sound is generated just prior to the arrival of the vortex at the 
step, within a distance equal approximately to the step height h. 


(3.8) 


Figure 3 illustrates how this picture is dramatically changed when account is taken of 
separation at the step. The incident vortex path is deflected sideways by the shed vortex 
(Figure 3a), and the two vortices proceed away from the wall in a direction initially inclined 
towards negative x. The times indicated on the trajectories correspond to those in Figure 2, 
i.e., t = 0 defines the instant at which the vortex T 0 would pass over the step in the absence 
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of shedding. The strength of the shed vortex increases as T 0 approaches the step until its 
value at Ut/h = l is about -1.2IV More details of the variation of T(t) are given below. It 
should be noted, however, that, because |T| > r o when Ut/h > 1, the vortex pair traverses 
a circular path which ultimately causes the vortices to return to the wall and subsequently 
to separate at a “reattachment” point far to the right of the step. In more realistic flows, in 
particular in the presence of mean flow over the step in the ar-direction, both the deflection 
distance away from the wall and the reattachment length would be much smaller. 

The sound pressure produced by this interaction is represented by the solid curve in 
Figure 3b. Also shown (dashed) are the acoustic pressures generated by the incident and 
shed vortices, which correspond respectively to the first and second terms in the square 
braces of equation (3.7). These separate contributions are both large, but of opposite sign, 
and interfere to produce a net radiated sound pressure that is much weaker than in the 
absence of shedding (shown dotted). 

When the initial stand-off distance of T 0 increases to d = h, Figure 4 shows that the 
trajectories and the pressure signatures are qualitatively similar to those of Figure 3; the 
acoustic amplitudes with and without vortex shedding are weaker, and the amplitude 
difference is smaller. If the vortex is assumed to be a crude model of a discrete eddy in a 
turbulent boundary layer, a larger value of d/h would be representative of turbulent flow 
over a small skin step, when the boundary layer thickness tends to be large compared to 
h. The case in which d/h = 4 is illustrated in Figure 5. Here the deflection of the incident 
vortex is relatively small (about three step heights), and the circular path of the vortex pair 
terminates just to the right of the step (near Ut/h « 3), the trajectory being reminiscent 
of reattachment profiles observed in experiments (Farabee & Casarella 1984, 1986, 1988): 
the “eddy” formed by the vortex pair is torn apart, T 0 proceeding to the right along the 
wall and away from the step, while T is captured within a “separation bubble” . Actually 
dT/dt changes sign at Ut/h ~ 5.2, at which time it is “released” from the edge S with a 
final circulation T = — 1.27IY The vortex then proceeds back towards the step and proper 
continuation of the calculation requires that a new vortex be released from S. However, this 
latter phase of the motion will not be pursued (being similar to the interactions discussed 
below for the backward facing step), since the interaction of the incident vortex T 0 with the 
step is now complete. The acoustic pressures with and without vortex shedding are seen in 
Figure 5b to be of very similar amplitude, but practically opposite in phase. 

The variation of r(t)/T 0 when d/h = 4 is shown in Figure 6. The vortex is released at 
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Ut/h « 5.2 as indicated. In all of these cases the sound is generated principally over an 
interval of time ~ h/U, with characteristic wavelength ~ h/M > h at the small Mach 
numbers for which the present theory is applicable. 

3.4 The backward facing step 

Typical interactions for a backward facing step are illustrated in Figures 7 and 8, 
respectively for d/h = 0.5 and 2. In the absence of shedding, the vortex paths and the 
acoustic pressure signatures are identical with their respective counterparts for the forward 
facing step when the time direction is reversed. This is because of our assumption that 
the span £ of the wall “wetted” by the unsteady flow is small compared to the acoustic 
wavelength. As £ -» oo, an observer at a given far field point x first receives sound from the 
interaction of the vortex and step occurring at the closest point of the step (corresponding 
essentially to the pressures calculated in this chapter), but this is subsequently augmented 
at progressively increasing values of the retarded time with sound generated by interactions 
at increasing spanwise distances on the step from the observer. 

When vortex shedding occurs at a backward facing step, the incident vortex is displaced 
temporarily away from the wall (cf. Figures 7a, 7b and 8a), thereby reducing the intensity 
of the generated sound. The incident and shed vortex form a vortex pair (with T « — T 0 
when Ut/h ~ 1) whose circular trajectory approaches the wall “downstream” of the step, 
where the pair split up as the flow “reattaches”. The incident vortex continues along a 
path parallel to the wall and the shed vortex translates back towards the step within the 
separation “bubble”. This vortex is “released” before reaching the step (respectively at 
Ut/h « 2.8 and 2.2 in Figures 7 and 8), whereupon a new vortex must be regarded as 
shed from S in order to continue the calculation. This secondary shedding has not been 
modeled because it occurs long after the interaction of r o with the step, and is probably 
not important in practice, since at these times surface interactions would be dominated by 
other dissipative mechanisms within the separation zone. 

When d/h = 0.5 the acoustic pressure (Figure 7) is much smaller than in the absence 
of shedding, and is “phase advanced” with respect to the no-shedding pressure signature. 
The separate contributions to the radiation from the incident and shed vortices are again 
both large, but of opposite sign, and their interference determines the amplitude of the 
sound. Figure 8 shows that (in contrast to the forward facing step) an increase in the 
initial stand-off distance to d — 2h leaves the relative amplitudes of the acoustic pressures 
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with and without shedding effectively unchanged until the shed vortex is released. The 
subsequent sharp increase in the radiation after the release of the shed vortex at Ut/hzz 2.2 
is of no practical significance, because it represents the sound generated when T returns to 
the step after being released. 
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4. CONCLUSION 

The unsteady drag exerted on a stationary body in incompressible flow can be expressed 
entirely in terms of the vorticity distribution. To calculate drag fluctuations produced by 
impinging turbulence it is important to include contributions from vorticity generated at 
the surface of the body. At low Mach numbers the drag is equal to the strength of the 
aeroacoustic dipole that dominates the acoustic radiation. The numerical results given in 
this chapter for an idealized vortex-step interaction suggest that vorticity in the separation 
zone significantly modifies the radiation compared to predictions based on a potential flow 
“edge scattering” modeling of the surface generated sound. This may be representative of 
the kind of interaction that occurs when a large scale, turbulent structure convects over 
forward or backward steps. 


d/h 

Attenuation (dB) 
forward step backward step 

0.5 

11.8 

10.1 

1 

6.3 

9.7 

2 

3.5 

9.0 

4 

1.4 

4.8 


Table 1. Attenuation of vortex-step interaction noise by vorticity production 

The influence of shedding from a step of height h depends on the stand-off distance d of 
the vortex from the wall prior to the interaction. Predictions for large or small values of 
d/h might be expected to be relevant to sound production by turbulent flow over a step 
when the corresponding ratio 6/h of the boundary layer thickness 5 to the step height is 
large or small. Vorticity production affects both the amplitude and the phase of the sound. 
An overall estimate of the attenuation of the step noise by the surface-generated vorticity 
can be made by comparing the respective net radiated acoustic energies £ s and E 0 , say, 
with and without vortex shedding. This is done by comparing corresponding values of 
S—oo* P 2 (x, t) di, where the interaction of the step with the incident vortex is regarded as 
being negligible beyond a “cut-off” time t max , which is introduced to eliminate spurious, 
large amplitude contributions from secondary interactions of the shed vorticity with the step 
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(that are suppressed in practice by viscous dissipation and additional vorticity production 
at the wall). This comparison is made in Table 1, where the attenuation -10 x log 10 (£ s /£ o ) 
(dB) is given for several values of d/h. The attenuation decreases with increasing d/h 
(“boundary layer thickness”), but less so for the backward facing step, indicating that 
“reattachment” on top of a forward facing step is noisier for a given value of d/h. 

For a forward facing step, when d/h is large, Figure 5 and Table 1 indicate that the 
amplitudes of the sound with and without shedding are roughly equal, but the phases are 
reversed; phase reversal is also evident in Figure 8 for the backward facing step. These 
results suggest that predictions of the frequency spectrum of rough wall boundary layer 
noise based on the potential flow scattering of the turbulence blocked pressure by very small 
roughness elements (see, e.g., Howe 1989b) are probably satisfactory for practical purposes, 
and justifies in addition recent estimates of the noise generated by turbulent flow over small 
fuselage skin steps (Howe 1997). 
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Figure 1 . Configuration of the incident and shed vortices for a forward 
facing step. 
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U[t]/h 


Figure 2. (a) Trajectory of incident vortex in the absence of shedding for d/h = 
(b) Non-dimensional acoustic pressure p(x, t)/[p 0 U 2 M(d/h) 3 (£/\x.\) cos#]. 


0.5. 
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U[t]/h 


Figure 3. (a) Vortex trajectories when d/h = 0.5. (b) , non-dimensional 

acoustic pressure p(x, t)/[p Q U 2 M(d/h) 3 (£/\x\) cos#]; •••••, acoustic 


pressure in the absence of shedding. The broken curves are the separate 
contributions from the incident and shed vortices. 
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U[t]/h 


Figure 4. (a) Vortex trajectories when d/h = 1. (b) , non-dimensional 

acoustic pressure p(x, t)/[p 0 U 2 M(d/h) 3 (£/\x\) cos#]; •••••, acoustic 


pressure in the absence of shedding. The broken curves are the separate 
contributions from the incident and shed vortices. 
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Figure 5. (a) Vortex trajectories when d/h = 4. (b) , non-dimensional 

acoustic pressure p(x, t)/[p 0 U 2 M(d/h) 3 (£/\x\) cos#]; •••••, acoustic pres- 
sure in the absence of shedding. 
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Ut/h 


Figure 6. Shed vortex strength r(t)/r o for a forward facing step when d/h = 4. The 
vortex is “released” from the edge S of the step at Ut/h ~ 5.2. 
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U[t]/h 


Figure 7. (a) Vortex trajectory for d/h = 0.5 with no shedding. 

(b) Trajectories with shedding. 

(c) , non-dimensional acoustic pressure 

p(x, t)/[p 0 U 2 M(d/h) 3 (i/\x\) cos#]; 

• ••••, acoustic pressure in the absence of shedding. 


108 




Report No. AM-97-027 


Boston University, College of Engineering 




U[t]/h 


Figure 8. (a) Vortex trajectories when djh — 2. (b) , non-dimensional 

acoustic pressure p(x, t)/[p 0 U 2 M(d/h) 3 (£/\x\) cos#]; •••••, acoustic pres- 
sure in the absence of shedding. 
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